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(57) Abstract 

A computer method manages risk in multiple-parameter physical 
systems perfomitng interrelated activities, where at least one of such 
activities is risk-related in that it may iiave ao outcome level which 
may fall outskle of boundary limits. The method establishes a course 
of action for the physical systems that facilitates preventing any outcome 
levels for risk--re]ated activities from falling outside of boundary limits 
(2S). The method assuntes the existence of a computational multiscenario 
decisioTMnaking model (17) that describes the physical .systems and 
deteimines, under some set of criteria, both feasible and desirable levels 
of their activities. The method finds a set of satisfiaUe boundary limits 
in computer nicmoiy (45), develops in computer memory a multitude of 
candidate strategies (75) that satisfy tlieae limits, describes the strategies 
in computer memoiy in fomuits of multidimcnsioinl outcome and regret 
matrices and jointly applies to such matrices multiple optimization criteria. 
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pOMPUTER AIDED RISK MANAGEMENT IN MULTIPLE- gARftMBTER 

PHYfiTHAL SYSTEMS 
Pield Of eh^ Invention 

This invention relates to caacputex aided management 
5 of risk in multiple -parameter physical systems having one 
or more risJc- related activities and more particularly, 
although in its broader aspects not exclusively, to 
conputer aided management of risks arising from decisions 
made under conditions of imcertainty. 

10 BackaiTQunrl nf the Invention 

A need for risk management decisions arises in a 
broad range of technological i industrial, and financial 
areas. Typical examples include operation of 
manufacturing, storage, and transportation facilities in 

15 an industrial logistics system, control of product mix at 
a factory, deployment of industrial equipment, electri<:al 
engineering and mechanical engineering problen^, 
inventory control, advertising campaign management in a 
marketing program, and management of a portfolio of 

20 financial assets, to name just a few. 

As an example, consider the problem of efficient 
lighting in a commercial facility. The facility ovmer 
has to provide the required lighting conditions in the 
building. He would like to Install state-of-the-art 

25 lighting systems and controls In order to curb energy 
consumption at this facility. The state-of-the-art 
equipment commands a premium price, mitigated however by 
rebate incentives from the local utility <:ompany. At the 
same time, the owner would like to minimize the cost of 

30 the lighting devices that are necessary to provide the 
required lighting conditions. 

In this example, the owner calculates ^he operatiiig 
costs of lighting the facility with dif4!erent types of 
devices, based upon data provided by the manufacturers. 

35 These data usually correspond to power factor at the 
facility being equal to 1.0- Real world conditions 
introduce a variability of power distribution <:o that 
facility that reduces the power factor and affects the 
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operation of Che selected equipment and, therefore, its 
operating costs. Such conditions induce a consideration 
of auxiliary devices that restore the power factor at the 
facility to unity. 
5 The risk management decision in the facility owner's 

problem is to find a combination of the quantity amd 
quality of lighting devices, controls, and auxiliary 
devices that minimizes the payback period but still 
affords protection from the twofold risks of both 

10 exceeding the planned operating costs and underachieving 
the desired lighting conditions. The decision should 
contain these risks within some acceptable limits. 

To demonstrate the universality of the need for risk 
management decisions, consider next a problem of 

15 efficient distribution of products by an industrial 

manufacturing company. The distribution system starts 
with the con5>any factories that manufacture the products 
and ends with buyers (such as wholesalers) who order the 
products. The system includes a network of distribution 

20 centers and warehouses, as well as transportation 

facilities to move the products. All of these may belong 
to the conipany, or may simply be used by it. All 
facilities of the system (namely, the factories, the 
distribution centers, the warehouses, and the 

25 transportation vehicles) are characterized by their 
production, throughput, or storage capacities. 
Similarly, iise of all these facilities Invokes their 
associated costs. If the distribution system Is not 
limited to a single country, costs and prices may need to 

30 be expressed in different currencies. Manufacturing, 

transportation, loading and unloading, and hauidling the 
products at warehouses - all of these procedures ^require 
resources and time. 

Demand depends on product prices which, in turn, may 

35 be related to cumulative product costs at the buyers* 
locations. It also depends on the behavior of 
competitors, which is determined partly by the con^ny's 
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own pricing and other policies , and partly on other 
factors, largely untaiown.' 

Any inability to meet the buyers' demand , and to do 
it on time^ involves explicit or inplicit economic 
5 penalties. Similarly » penalties arise if the ^conpany 

procures work force, equipment, and materials to meet its 
planned production targets, and then has to change its 
plans, causing mismatches. 

The values of all parameters of the distribution 

10 system, including its technologies, the needed production 
and transportation time, the required resources and 
capacities, demand, prices and costs, currency -conversion 
rates, and penaaties (especially inplicit penalties), are 
not known exactly. They may also change over time. The 

15 values for these parameters can only be estimated or 
forecast. The risk management part of the effi-cient 
distribution problem is to find a combination of 
technologies, production targets « inventozry levels, and 
transportation flow at all stages of the distribution 

20 system during the planning period, as well as of product 
selling prices and levels of demand to be satisfied, so 
that no production, warehousing, or transportation 
capacity constraints are exceeded and total profits are 
maximized, while the risks of insufficient profits or 

25 losses, penalties, foreign exchange rate changes, or 
\inmet demand and broken schedules are kept within 
acceptable limits. 

Finally, let us consider a financial portfolio 
mainagement problem. For simplicity, assume that the 

30 portfolio may include only fixed income securities of 
different maturities but of one general type, such as 
bonds issued by the United States Treasury. The 
portfolio does not include corporate and municipal bonds, 
stocks, financial instruments in currencies other than U. 

35 dollars, mortgage- based securities, or derivative 

financial instruments, such as options. 

In this last example, the portfolio manager has 
exact data about the conposition of his portfolio, that 
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is, about the face value of portfolio bonds, by issues. 
The manager also knows all characteristics of each 
existing Treasury bond issue, both present and not 
present in the current portfolio* These characteristics 
5 include the issuers date of maturity, <:oupon or discoimt 
rate, the schedule of coupon payments, the coupon 
interest that has accrued on the issue from the time of 
the laBt coupon payment, the transaction costs on 
acquiring or selling the bond, callability, and the 
10 availability of the issue for purchase. 

In this example, the manager also knows ^he latest 
bond market quotes on bid and ask prices for all existing 
Treasury issues, although these quotes may differ from 
the real execution prices of bond trading. This 
15 difference may exist even if the bond is traded 

(purchased or sold) at this very moment, especially for 
bonds not actively traded in the market- However, both 
the bid and ask quotes and the execution pri<:es for each 
bozxl issue depend upon the supply/demand relationship for 
20 that issue, which changes all the time. Therefore, if a 
bond is purchased or sold not immediately but later, this 
relationship may change drasti*cally, entailing the 
corresponding price changes. For any time in the future, 
the portfolio manager does not know in what direction and 
25 how much prices will change for any issue. Moreover, 
changing bond prices affect not only new trades: the 
worth of the whole portfolio is regularly re-evalua-ted 
(marked* to-market) at current pri<^es. 

In this exaii5)le, the portfolio manager wants to 
30 msocimize portfolio returns. However, he also has to 
carry out certain obligations ^o the portfolio owners 
(Investors) . Perhaps the most in^rtant obligation is to 
make scheduled payments to investors - either some 
contractually specified amounts, amo\ints that stand for 
35 the returns on investment and repayment of the investment 
principal, or amounts that symbolically represent 
advances on the Investment returns that are expected in 
the future. 
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Besides these payments « the portfolio will have in 
the future some other cash inflows and outflows. The 
inflows are mostly new investments in the portfolio* 
coupon payments from the Treasury on the portfolio bonds, 
5 and the principals of the portfolio bonds that have 

matured. The main outflows are the withdrawals by the 
investors from the portfolio and the administrative costs 
of portfolio management. Wiile the coupon payments fully 
depend on the conq^osition of the portfolio and the 

10 management costs can be anticipated with sufficient 
accuracy, both the new portfolio investments and 
portfolio withdrawals can at best be "educated guesses •** 
The portfolio can trade bonds, i.e., sell the bonds 
currently in the portfolio and, using these funds 

15 together with new inflows, purchase some other bonds. 

The permitted portfolio trading activities are 
restricted by a number of laws, rules and constraints of 
fiduciary, regulatory, tax, and other origin. There may 
be constraints on borrowing, margin -trading and other 

20 leveraging of the portfolio funds, short sales, and so 
on. One of the main constraints is a fiduciary 
requirement that commonly obliges the manager to preserve 
the principal capital of the portfolio, that is, to 
protect the portfolio against unacceptable risks. 

25 The portfolio manager can base his portfolio 

decisions either solely on the latest bond market quotes 
which he knows, or on a combination of these quotes and 
future bond prices that can be expected for some moment 
of time yet to come. 

30 The risk management problem of the portfolio manager 

is to find a planned ^onibination of bond trades so that 
all constraints on the portfolio activities are met and 
the portfolio returns are maximized, while the portfolio 
funds are protected against losses ^hat would exceed the 

25 acceptable risk limits. 

These three examples demonstrate the extreme 
complexity of making risk management decisions. Still, 
these are relatively simple situations - -real life 
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decision-xnalclng in financial and industrial business 
organizations can be much more con^licated. For 
instance r the financial portfolio in ^he last exanqple 
could include not only Treasury bonds, but also -corporate 
5 and municipal bonds, stocks « futures « options, finstncial 
instruments in other currencies, and other ^ypes of 
securities. 

It should be noted that the risk roanag^nent problems 
discussed above are real physical problems arising in 

10 real physical systems. (A portfolio of fixed ixicame 

securities is also a physical system - a sec of physical 
bonds.) While it is true that this invention represents 
significant quantitative aspects of these phy8i<:al 
problems by mathematical models, the purpose of th^se 

15 models is to make decisions about target values which are 
then used in the physixtal world ^o construct or operate 
physical systems. Typical prior art -examples of such 
mathematical models for making decisions about physical 
systems are the use of linear programming <LP) for 

20 efficient resource allocation or for optimizing system 

operational parameters, the use of scenario optimization 
for the management of a portfolio of financial options, 
and the use of mathematical equations to construct radio 
antennas or to control riibber- molding operations • 

25 Moreover, the mathematical models used in this 

Invention are, as a rule, too coxnpllcatied for appl ideation 
in reasonable time without a computer. Therefore, their 
use involves changing the physitral condition of the 
conputer memory, thus virtually creating a new state of 

30 the computer. 

To exercise risk management, it is first necessary 
to define "risk." 

Every decision in any of the areas listed above, 
such as electrical engineering, industrial logistics, and 

35 portfolio management, involves, on the one hand, a 

specific allocation of resources and, on the other hand, 
specific outcomes from different activities of the 
physical system. These results may be as diverse as 



wo 98^3776 



7 



PCT/USy7/16446 



coats or profits « returns on financial portfolios, 
quality of products « the composition of product mix* 
consiunption of electricity, penalties for brolcen 
schedules or imsatisf ied demand, or the amount and 
5 chemical conposition of waste water* The values of some 
of these results are of special concern to decision- 
makers (DMs) , who consider such results undesirable or 
even potentially dangerous and want them either to scay 
within B<mB boundary limi-ts (which may or may not be 

10 known in advance) , or to be as low or as high as 

possible. The systems activities that are subject to 
special concern will be called "risk-related activities*" 

Outcomes in risk*- related activities depend both on 
the targets to be achieved and on the allocation of 

15 resources within the physical system. In general, it is 
this dependency that is described by a mathematical 
model. As a rule« the parameter values of the dependency 
relationship in the system, or of the mathematical model, 
are not known with certainty and ^can, at best, be 

20 estimated or forecast, and often just guessed. Even if 
they are forecast, it is still hardly ever, if at all, 
possible to obtain a reliable forecast, especially for 
the long tenii. The consequences of a decision about the 
allocation of resources in a physical system are, 

25 therefore, uncertain and involve ^he possibility of 

unesqpected and, possibly, undesirable or even dangerous 
outcomes • 

"Risk" is defined here as a magnitude of outcome 
levels of undesirable or potentially dangerous activities 

30 that have fallen outside the relevant boundary limits. 

Risk management, then, is the capability to estimate, to 
avoid, to control or to reduce the extent of such 
occasions and reduce the probability of their 
occurrence* "Risk" can be defined in many ways, both 

35 absolute and relative, and this invention can be applied 
in the framework of any of these definitions* The exact 
definition of risk is here irrelevant, so, for 
simplicity, only one definition is used. 
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To understand the need for the present invention, 
consider the current state-of-the-art both in the general 
area of decision-making and, specifically, in risk 
management * 

5 Roughly, from the time the computer era began - that 

is, in the late 19408 and early 19508 - great advances 
were made in two fields important to decision-making: 
Operations Research/Management Science (OR/MS) and 
Decision Science (DS) • OR/MS primarily deals with 

10 optimization models, while DS analyzes alternative 

strategies under uncertainty and selects one candidate 
strategy over all others* 

A rational, natural and customary way to make a 
decision is to recognize the uncertainty of the future 

15 and the lack of knowledge about the present, to represent 
the uncertainty and lack of knowledge through scenarios, 
and to consider the outcomes - possibly in many 
activities - under each scenario, given a course of 
action or strategy. Arranging information about these 

20 outcomes into an "outcome matrix** is a good technique tor 
systematic analysis of the data. This is the basic 
approach of cla88i<:al DS, which starts from an out<:ome 
matrix, or, more specifically, often from a "payoff 
matrix" - a special case where all outcomes are 

25 quantified and are of the same type, such as profit. 

However, this assumes that scenarios, candi<3ate 
strategies, and "strategy versus scenario*" out<:oine8 are 
specified beforehand. Thus, classical DS in effect is 
"passive" - it withdraws itself from tasks that 

30 constitute 95 to 99 percent,. or even more, of the total 
effort* Instead of addressing the whole real world 
problem of making a decision, DS limits itself to the 
last, and often the easiest, part of the process. 
In contrast, OR/MS uses "active" optimization 

35 models. These models address the main part of the 

decision-making problem not covered by DS: they formulate 
a plan or a strategy. Using mathematical equations and 
inequalities, t:he models define a region of feasible 
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solutions (i.e., solutions that do not violate those 
equations and inequalities) and find in that region a 
solution that is "the best" from the point of view of one 
or more criteria of optimality. 
5 The most widely used optimization model is an LP 

model. It is also the basis for and a major -component of 
all other more sophisticated mathematical programming 
models « such as integer or non- linear programming models. 
It defines a set of interrelated activities and a set of 

10 constraints on the level of each activity, as well as on 
the levels of seme specified linear functions of these 
activities. The LP model also defines the "objective 
function" as the total sum of the net benefits ^benefits 
minus expenses) of activities, which also is a linear 

15 function of activity levels. The solution of an LP model 
is based on a single criterion: finding the allocation of 
resources to maximize the value of the objective 
function. One method of solving an LP model is described 
by United States Patent No. 4,744,026 to Robert J. 

20 Vanderbei, United States Patent No. 4,744,027 to David A. 
Bayer et al., and United States Patent No. 4,744,028 to 
Narendra K. Karmarkar, all issued May 10, 198B. 

When applied correctly, LP models have many valuable 
advcuitages. The models can integrate, connect, 

25 coordinate, balance, and jointly analyze different 

factors, operations, and territorial or functional parts 
of a physical system. They can find hidden opportunities 
for improvement and are easy to set up, although this 
simplicity can sometimes be very deceptive. They ^an 

30 also derive plans from the initial data, without the need 
of losing time and effort on intermediate -analyses, and 
thus increase the speed of decision -making. 

LP models became invaluable tools for dealing with 
"closed" decision-making problems, problems that exclude 

35 any significant deviations from the status quo in 

important decisions. In the petroleum industry, whioh 
began to apply LP in the late 1940s, mathemati^cal 
programming concepts penetrated all facets of short -term 
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busiiieas planning, from supply, distribution and refinery 
planning to product blending and process -control. 
Scheduling of Air Force planes was a triuncphal 
application of sophisticated, large-scale integer LP 
5 models during the 1991 Gulf War. 

Conventional optimization techniqxies have never been 
successful, however, in dealing with "open" problems that 
predominate in long -term and strategic planning. This is 
because prudence and moderation, the -two crucial 

10 components of mature decision- making and risk management, 
are not among the advantages of optimization models. LP 
models seek extremes and are stopped only by such model 
constraints as equations or inequalities. Bvaen a 
minuscule alteration of input data, well within the 

15 meucgin of possible error, may cause a change of solution. 
Moreover, the solution always switches from one extreme 
to another, so that the resulting change <:an be 
disproportionately large. Therefore, solutions of LP 
models are inherently xrns table; they introduce an 

20 additional risk component of their own. 

InsteQDility of solutions causes three major 
difficulties in applying such models. The first such 
difficulty comes from uncertainty. A model may provide 
valid results only if the model data are sufficiently 

25 accurate (which means the modeler must have adequate 

knowledge of both the present and the future) , or if ^he 
major decisions in the optimal solution remain 
sufficiently stable as data varies regarding the possible 
actions and their consequences. Second, LP models have a 

30 simplistic, well defined, one-dimensional goal, while the 
nte* goals usually are more diverse, conflicting, and 
ambiguous. Third, by definition, the models are 
incomplete and they may omit in^rtant factors, 
considerations, and constraints, such as long-^erm 

35 considerations in a short-term model. 

In other words, LP models fully confirm an 
observation of Oscar Morgenstern ("On the Accuracy of 
Economic Observations," Princeton University Press, 1950, 
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p. 45) that every type of numerical observation, 

based upon a mathematically formulated model, imposes 
restrictions upon the data. If these restrictions cannot 
be met, the operations become iopossible, even if the 
5 underlying model should be free from objections." 

Remedies have been proposed to deal with some of 
these drawbacks. For example, stochastic programming 
(SP) , an inqportant mathematical programming extension of 
LP. is intended to deal with uncertainty. To a certain 

10 degree, SP performs this function, but it of^f^ers only an 
ixnplicit and limited protection against rislc attendant to 
uncertainty. Its risk protection is valid only in a 
statistical sense: if ita optimal solution is repeatedly 
implemented a very large nuinber of times, it will 

15 eventually prove the best. However, with a few 

exceptions, the decision-making situations are either 
unique and non- repetitive or are repeated only a small 
number of times. If a harmful or adverse situation 
happens during one of those times, the losses resulting 

20 from the optimal" decision may never be recouped. 

Also important, SP does not offer choi<:e8 to <:he 
deci8ion*makers: it constructs a single solution and 
declares it to be the optimum. SP therefore deprives 
them of enormous advantages coming from the use of 

25 outcome and regret matrices and DS criteria. 

However, the crucial point is that SP usually is 
simply inappropriate as a basis for making complex 
decisions in the relevant fields. In "Risk, Uncertainty 
and Profit," by F. H. Knight, University of Chicago 

30 Press, 1921, a clear distinction is made between 

"insurable risk" and "non- insurable uncertainty." In 
that approach, insurad>le risk is said to exist when the 
probabilities of outcomes are known exactly and aire 
derived on an objective basis; that is, they are 

35 calculable on the basis of relative frequencies or 

similar data. Non- insurable uncertainty exists in the 
absence of objective and known probabilities. 
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At the same time, probabilities are considered to be 
calculable and adequate only for "repetitive phenomena of 
a standardized variety such as occur in games of chance, 
in actuarial science, in genetics, and in statist i<;al 
5 mechanics." ("Decision Analysis - Introductory Iiectures 
on Choice under Uncertainty" by HoiRrard Raiffa, Addison 
Wesley, 1970, p. 274). In contrast, the decision-making 
problems to be addressed by SP deal with complex 
economic, financial, technical and social phenomena, 

10 which, as indicated above, are non- repetitive. At best, 
the probabilities related to these phenomena include 
subj^tive judgments and more or less educated guesses. 
In these problems, probabilities are t:he Least reliable 
part of the input data. 

15 Recognizing this need in probabilities of future 

events, a leading practitioner of SP states, for 
instance, in the united Stages Patent No. 5,148,365 <o 
Ron S. Dembo, issued September 15, 1992, that "For those 
of skill in the art of portfolio management, the 

20 probability of the various scenarios can be guesstiroated 
with reasonable accuracy based on experience" (column 8, 
lines 56-59) . Notably, the author does not even mention 
objective probabilities which, as indicated above, are 
the precondition for proper applications of SP. He would 

25 be quite satisfied with the "guesstimated" probabilities, 
subjectively assumed by the portfolio managers. 

However, even these lo%*er'ed data demands cannot be 
met. Those "of skill in the art of portfolio management" 
disagree with the author's high evaluation of their 

30 capabilities and consider it to be little but wishful 

thinking. A leading financial forecaster writes: "To be 
sure, most forecasters* expectations do not work out at 
all- For instance, in June 1990, 88% of economists 
predicted continued economic expansion for at least a 
35 year. A month later, the worst recession in a decade 

began. As merely the latest example, a June 1994 survey 
of 29 of the coxmtry's most influential money managers 
showed that ^11 of them expected the long bond yield to 
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remain below 8% during the rest of last year, it was 
above 8% three months later. Evidence of the failure of 
conventional forecasting methods is more than anecdotal. 
According to The Wall street Journal, a study of its own 
S surveys since 1982 of the country's top economists 

reveals that in the aggregate, these acknowledged experts 
predicted accurately the direetion (forget the extent) of 
interest rates only 2S% of the time, which is half the 
success rate that would be produced purely by guessing." 

10 (Robert R- Prechter, "At the Crest of the Tidal Wave« " 
New Classics Library, 1995, pp. 19-20) • The latest 
surveys provide similar results. fThe Wall Street 
OfiUmal. August 6, 1996, p. AIS) . 

Moreover, if the applications of SP are 

15 controversial even when they Involve just the insurable 
ris)c, because the "optimal" strategies do not 
sufficiently protect from risk, then, under uncertainty, 
when the objective probabilities are not known, this 
method becomes even more controversial. Therefore the SP 

20 procedure has to be based on a combination of two 

inqplicit premises. The first aseiunes t:hat there exists 
an objective optimum under uncertainty, that it can be 
found by objective methods, and that it is just a matter 
of technique to find that optimum - namely, a matter of 

25 applying the correct model and getting correct data. The 
first premise also assumes that these techniques have 
been sufficiently attained in modeling and solving the 
problem under consideration.. The second premise concedes 
that the optimum may indeed be subject to qualifications, 

30 such as the personal risk attitude and subjective 

preferences of the dec is ion -maker, but assumes that these 
still can be incorporated into the SP model on the basis, 
say, of prior observations of the DMs* behavior and 
attitude, and that, again, these requirements are met in 

35 the problem under consideration. 

The first assumption has been proven wrong. Both 
the authoritative "Games and Incisions** by R. Duncan Luce 
and Howard Raiffa, John Wiley & Soils, 1957, pp« 274-303 
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and 324-326, and later literature, euch ae "Decision 
Malcing under Risk and Uncertainty: Mew Models and 
Enpirical Findings" ed. by John Geweke, Rluver Academic 
Publishers, Dordrecht, the Netherlands, 1992, pp. 1-10, 

5 show that, even in the sinplest case of a two-dimensional 
payoff matrix and a single decision outcome to be 
considered, such as profit, there are several DS criteria 
for decision-making under uncertainty - that is, methods 
for comparing and selecting strategies. (For instauice, 

10 DMB can base their choice on the average profit, the best 
case profit, the worst case profit, and some combination 
of the above . ) None of the known or even conceivable 
criteria of DS is perfect or "the best*" Bach has 
faults, such as violations of transitivity, that are 

15 revealed under scsne specific conditions. Thus, even in 
the slxr^lest case, it is impossible to make the best 
decision under \incertainty in a general, unique, 
objective, and theoretically correct manner. 

The second assumption is unrealistic, at least in 

20 the foreseeaUt>le future. As shown, for exanple, in -the 

above cited "Decision Making under Risk and Uncertainty: 
New Models and Empirical Findings" pp. 11-16, the 
existing theories of personal choice under uncertainty, 
such as the e3qpected utility theory, are still evolving. 

25 They cannot yet deal successfully even with some quite 
sinple but paradoxical decision -making situations. If 
and when this process successfully ^ends and some 
comprehensive and consistent decision theory, both 
normative and descriptive, becomes a reality, its 

30 "attitude extracting" procedures still are likely to be 
lengthy, cumbersome, imprecise, and impractical, not 
suitable for real life decision-making, especially in 
complex business situations that require reasonably quick 
decisions. (The expected utility theory lays claims only 

35 to normative correctness but not to a good descriptive 
characterization of choice \xnder risk and uncertainty. 
Therefore, it is doubtfxil that the theory can provide a 
satisfactory "attitude extracting.") 
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The **robu8t optimization" (RO) approach described* 
for instance, in "Robust dptimization of Large- Scale 
Systems** by J« M. Mulvey, R. J Vanderbei and S. A. 
Zenios, Operations Research, v. 43 <1995), No. 2, pp. 
5 264-281, is an extension of SP and has some advantages 
over SP: it allows solutions that are relatively stable 
(that is, the optimal solutions under different scenarios 
remain closer to each other than under SP) • Also, RO is 
multicriterial and allows tradeoffs between several 

10 criteria of optimality. However, similarly -to SP, RO 

still relies on scenario probabilities being objectively 
known. (In the quoted article, the authors -circumvent 
this crucial issue by siiqply mentioning, in passing, on 
p. 265 "the probability of ^he scenario." They do not 

IS explain how they succeeded in obtaining these 

probabilities and whether these are objective or not.) 
Moreover, even If these unrealistic expectations are met, 
RO would again offer only long run, "statistical" 
protection from risk, while its short-term results may be 

20 disastrous. Although, by changing the weights of 

different criteria, RO can form several strategies, it 
neither constructs payoff or outcome matrices nor applies 
them for conparlson and selection of the best strategy. 
Also, RO does not use clustering and therefore has to 

25 solve models with enormous nxambers of scenarios . 

Finally, too much importance is attached to meeting the 
initial constraints of the model (see later) • 

All above considerations about OR/MS and its tools 
refer to the first part of the decision -making process, 

30 namely, to the formation of candidate strategies. As for 
the second part (selection of a strategy) , this is the 
province of DS. It was indicated above that there are 
several DS criteria for decision-making under 
uncertainty. Most of th&x\ are based on the "strategy 

35 versus scenario" payoff values for a stratiegy - either on 
individual values, such as the best: payoff and the worst 
payoff of a strategy, or on values derived fr-oro 
individual payoffs, such as the average payoff. Three 
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best« most con^rehenslve and sophletlcaeed ^ri^eria of 
payoff type are the optimisni-pessiinism index (OP) 
criterion, the partial ignorance (PI) criterion, and the 
modified insufficient reason (IR) criterion. These three 
5 criteria, which are previously ]cnown, provide the basis 
for new methods of this invention. 

All three are "synthetic" criteria, which means that 
they are quite general and include as special cases 
other, simpler criteria. For instance, both minimax and 

10 maximin payoffs are special cases of all three criteria, 
the expected payoff is a special case of the PI 
criterion, etc. The synthetic criteria perform an 
extremely important role: they minimize ^he negative 
impact of absence or lac)c of knowledge about 

15 probabilities of future events and their <:ombinations 
(scenarios) . 

As mentioned before, probabilities are the least 
reliable part of input data, and decision analysis under 
risk and uncertainty cannot therefore generally dispense 

20 with subjective judgments, including judgments on 

probabilities. The goal is to minimize both the impact 
of these judgments and the effort required to form them. 
It is especially iB5)orcauat to arrange the introduction 
and use of probabilities and other judgments in a manner 

25 least detrimental to successful decision -making, which 
means to postpone their use until the latest possible 
stage of analysis. 

Fortunately, probabilities do not have to be used at 
the initial stages of the analysis, as is done in SP, 

30 decision tree methods, and so on. They also <can be 

compressed into a very few parameters that estimate -the 
overall degree of uncertainty. Furthermore, data 
requirements can be made less stringent by allowing the 
values of these parameters to fall within broad 

35 intervals, rather than correspond to a single value. 
Synthetic criteria meet all these requirements. 

All new methods are also synthetic and tiherefore 
have the same advantages. However, these .methods 
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additionally combine OP, PI and IR wit:fa the concept of 
"regret" introduced into Decision Science by Savage in 
"The theory o£ statistical decision, " Journal of the 
American S t^atiatical AsBOCiation, 46 (1951), pp. 55-67. 
5 Regret is basically a cost of uncertainty; it is derived 
from the payoff matrix and characterizes the risk, or 
regret, or opportunity lost because of choosing a wrong 
strategy. Regret may also be defined as the potential 
for reducing risk by switching strategies. See 

10 calculation of regret by Bq. (7) . In some special 

connotations, but not in this invention, "regret" is 
defined as the difference between a given benchmark, such 
as the performance of the stock market, and the actual ^r 
projected results. In previous state-of-the-art Decision 

15 Science methods, regret has been used only in sinple non- 
synthetic criteria, such as average regret or minimax 
regret (see, for instance, J. R. Buck, "Economic Risk 
Decisions in Engineering and Management," Iowa State 
University Press, 1989, pp. 313-334) . To the best of my 

20 knowledge, the con?>rehensive OP, PI and iR criteria have 
not been applied to regret; it is done for the first -time 
in this invention. 

In addition, the new regret-based synthetic methods 
naturally give rise to decision formulas and graphs that 

25 use and expand the concept of "efficient frontier" {see 
later) . The proposed new methods are invaluable €or 
finding desirable limits on tightening the discretionary 
constraints. 

A crucial dif f i<nilty in applying DS is that none of 
30 strategy selection criteria is "the best" under all 

circumstances . This opens the way to ^combining criteria 
* another method of this invention. 

The need to deal jointly with the totality of 
coxnpli-cated decision -making iss\ies such as uncertainty. 
35 the multiple criteria involved in real world decision- 
making, and the incompleteness of mathematical models, 
creates additional difficulties. Especially crucial is 
inconpleteness : it cannot be eliminated in principle. 
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since the only complete model of a reality Is t:he reality 
itself. 

This analysis shows that both general approaohes ^o 
decision-maKing (OR/MS and DS) have serious flaws if used 
5 as mutually exclusive tools, as is the current practice « 
Let us see how these general considerations are reflected 
in the specific field of risk management. 

Both the theory and practice of risk management are 
most advanced in the financial industry, therefore we 

10 begin our survey there, although the issues and 

techniques that are specific to managing portfolio risks 
are not directly addressed in this patent application. 

Modem Portfolio Theory constructs the efficient 
frontier - the risk/ return curve, whi^h defines a 

15 portfolio with the highest expected return for a given 
level of risk, or the lowest level of risk for a given 
level of expected return. In other words, it attempts to 
"optimize" a portfolio, to find the best tradeoff between 
expected return and expected risk. Iiftiat is "^he best** is 

20 determined by the subjective risk attitude of ^he 

decision-maker- However, in practice <his approach is 
"passive": it does not form portfolios but ranks only 
"external" portfolios that are developed outside the 
system. Moreover, the «9eakest link of this approach is 

25 its inability to evaluate and manage risk sufficiently 
well. Depending on market conditions, the for«:ast 
levels of risk may prove to be good or bad approximations 
of reality. 

Until the early 1980s, risk management in financial 
30 institutions was mainly limited to the use of 

Asset/Liability Models (ALM) • That methodology estimates 
future earnings under a number of probable scenarios of 
economic and financial conditions, projects future cash 
flows for one or more candidate investment strategies, 
35 derives final assets and liaUoilities resulting from each 
strategy under each scenario, and presents the estimated 
returns for all "strategy versus scenario" combinations 
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In t,he format of a payoff matrix. ALM thus follovs the 
approach of DS« but its iinplementation has serious flaws. 

MM wae primarily intended for such institutions as 
commercial banks « where both assets and liabilities were 
5 relatively illiquid and were priced on an accrual basis. 
ALM is based, homver, upon a false assunption that gains 
or losses occur when they accrue. To find out returns, 
AIM needs simulation over long periods « until most 
portfolio transactions mature. Trading items, which must 

10 be marked- to -market, are treated separately, and it is 

difficult to arrange hedging between trading and accrual 
items. For that purpose, "proxy values" (that is, 
approximations to market values) have to be estimated for 
accrual items. ALM also has other faults. Again, one of 

15 these is that ALM is passive - it provides no means for 
devising a good investment strategy and evaluates only 
"external" candidate portfolios. A second fault is that 
ALM is not capable of dealing with a large number of 
scenarios that might be needed because long -term 

20 scenarios are not accurate. The acciumilation of all 
these faults is worrisome. 

Unfortunately, during the past two decades several 
trends have evolved that made risk mauiagement in finance 
both more difficult and more necessary. Some of the most 

25 inqportant of these trends are: 

(a) Securitization of financial instruments, 
increase of their liquidity, wide use of more volatile 
instruments, such as derivatives and especially options, 
and moving from accrual accounting to frequent 

30 revaluation and marking-to-market of positions; 

(b) Increased volatility of financial markets, which 
is due to their globalization, advances in information 
technology, and growth of mutual funds, especially those 
specializing in emerging markets; 

35 (c) Increased trading, and especially the 

Institutional trading for an institution's own account; 
and 
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(d) Bntphasis on performance, which as a rule can be 
improved only by aseiiming higher risks, that sometimes 
lead to rogue trading, fraud, and eventual financial 
disasters . 

5 These trends cause concern about risk control among 

both the institutional managers and regulators. Two risk 
evaluation methodologies have been developed to meet 
these concerns. Value-at-Risk (VAR) considers risk that 
arises from random market movements, while Stress Testing 

10 deals with risk of the worst-case scenarios • Both 

methods have their advantages and disadvantages, amd they 
are best applicable under different circumstances. 

The VAR method assumes that rate and price movements 
of financial instruments can be described in a 

15 statistical fashion. If VAR is applied at times when 
this assumption is correct and the markets are 
statistically stable, the method provides an estimate of 
the loss that is expected to occur no more than, say, 5 
percent of the time. 

20 However, VAR depends heavily on estimates of 

volatilities and correlations that are derived either 
from historical data or from the values "in^lied" in 
current market prices. The trouble is that when a market 
collapses or makes a sharp move, that is, when we really 

25 need the risk control method -to work, all these estimates 
become irrelevant, because actual volatilities greatly 
exceed the estimated values • by at least several 
multiples . 

The Stress Testing method uses defined scenarios, 
30 including those for unstable markets . The scenarios <:an 
be simulated on the basis of both market -conditions for 
selected periods in the past and "educated guesses." 
This method provides more information of the expected 
portfolio performance, but it is computationally 
35 demanding even for a sharply restricted number of 
scenarios - 

These two methods can be used in combination, 
benefiting from the advantages of each. Their joint use 



wo 98/13776 



21 



FCTAJS97/I6446 



does not overcome, however, their coznmon fault: both 
methods ar *<paBsive" in that they do not necessarily 
generate good portfolios hut rather only evalixate 
portfolios constructed elsewhere* Both methods, and 
5 especially their combination, xian be used as a prelude to 
"optimising" portfolios, as defined above. 

There are attenqpts to provide a combination of 
portfolio optimization with risk management by SP. The 
trouble with these atteirpts is that SP protects from risk 

10 (and only from ■insurable risk") only in a stable market, 
where historical statistical parameters are valid* 

To sum up, no methods used in portfolio management 
are quite satisfactory or provide good risk msmagement. 
Outside the finance industry, both the theory and 

15 practice of risk management are much less advanced. A 

number of large companies follow the approach of scenario 
planning, which can serve as a foundation to risk 
management- As currently used, however, this methodology 
as a rule is wrongly focused on the definition of 

20 scenarios, rather than on rigorous develofxnent of 
candidate strategies, and especially compromise 
strategies . 

At the same time, the business world has become more 
volatile* Uncertainty has become the rule xather than 

25 exception, and it too often brings unpleasant suirprises. 
As in finance, this makes risk meuiagement both more 
difficult and more necessary. 

A common factor in all current practice is that, 
when dealing with multiple -parameter physical systems, 

30 effective risk management escalates quickly in complexity 
until it is literally beyond the capacity of the human 
mind to handle on any basis other than that of an 
educated guess or a "rule of thumb." The existing 
computational approaches are also inadequate. Moreover^ 

35 many of them, such as LP models, increase risk by adding 
a risk element of their own. There is, therefore, a real 
and continuing need for tools that will aid in valid risk 
mainagement . 
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SaBPftrv Qg the Inventipn 

This Invention deals vith decisions about the future 
activities of physical systexns that have potentially 
dangerous, risk-related activities. "Risk" is defined 
5 here as a magnitude of outcome levels of these activities 
that have fallen outside some boxmdary limits, which may 
be either predetermined, or determined in the process of 
making the decision, or both. This invention provides 
new techniques in five broad directions. 
IQ First, whatever model (used herein in -the broad 

sense of encompassing a formula or a method as well) is 
used to make the decisions, as a rule it is not always 
completely reliable for selecting the best decision. 
Such a model may, however, be used for an easier task - 
15 screening out bad and risky decisions. For that purpose, 
the model has to be complemented by adding e3q)licit risk- 
limiting constraints based on specific boundary limits. 
The "risk-limiting constraints" are defined here as 
either mathematical equations and inequalities that, once 
20 added to the model, eliminate or reduce outcomes falling 
outside the boundary limits, or other means (such as high 
fines and penalties) that play the same role, of limiting 
the risky outcomes. 

Second, the risk- limiting constraints in this 
25 invention are of two different types - predetermined 

constraints and discretionary constraints . Predetermined 
constraints are based upon generally recognized saCety 
and regulatory considerations, notions about the 
tolerable level of the enterprise's profitability, and so 
30 on. After these constraints .are met, this invention 
introduces discretionary constraints during 4:he 
conputational process of strategy formation - to find, 
analyze and fine-tune tradeoffs and to construct 
alternative candidate strategies. As a rule, 
35 predetermined constraints refer to the whole model, but . 
discretionary constraints are usually imposed on outcome 
levels of specific risk-related activi<:iee under specific 
scenarios . While the predetermined constraints may be 
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used In other models, discretionary constraints are 
unique to this invention. 

Thirds as a rule, the addition of risk- limiting 
constraints makes the expanded model more complex and 
5 thus necessitates the use of con^uters in finding the 
solutions of the eaqpsuoded models in a reasonable time. 
State-of-the-art optimization methods currently available 
for solving these complicated models, when applied 
without due attention to risk, are insufficient and often 

10 counterproductive, because any such method tends to add a 
substantial risk coitponent of its own. This invention 
provides new ways of dealing with that complexity in 
constructing candidate strategies (a "strategy" is a set 
or a sequence of possible decisions). These new ways 

15 include the use of multiple single- scenario optimization 
models, as well as both single and multiple multiscenario 
models. 

Fourths the special structure of <he multiscenario 
models, coxtibined with the use of a relatively small 

20 number of strategic variables and with the understanding 
that a good solution roust be good for worst scenarios, 
leads to the development through this invention of a 
novel, highly efficient and quickly converging 
decomposition method for solving large mathemati<:al 

25 programming models. 

Fifth, decision-making usually comprises two stages: 
construction of candidate strategies and selection of one 
of them, to be inplemented in the physical system or 
systems. This invention provides techniques to deal 

30 jointly with the whole decision- making process - not only 
with the first stage, but also with the second stage. At 
the second stage, these techniques provide new methods of 
comparing and selecting strategies. They also deal with 
group decisions, multicriterion decisions, qualitative 

35 attributes, and outcomes in multiple risk- related 
activities . 

Broadly, the invention takes the font o£ a con^uter 
method for managing risk in multiple parameter physical 
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systems performing lnterrelal:ed activities, where at 
least one activity is risk- related in that it has outcome 
levels that may fall outside boundary limits. Such 
boundary limits may be predetermined. The method 
5 specificeaiy helps manage risk by establishing a course 
of action (strategy) for the physi<:al systems preventing 
any outcome levels for the risk related activities from 
falling outside these boundary limits. The method 
assumes the existence of both boundary limits sucxd a 

10 decision-making model chat defines, under some set of 
criteria, the best solution of the model and the 
desirable levels of the activities of the physical 
systems in that solution. The initial model may be 
stored in coaqputer memory either directly in its 

15 computational form, or in its symbolic, mathematical 
form. In the latter instance, the model is converted 
into a computational model in computer memory by 
replacing symbols with nxmibers (i.e., the values of 
parameters) • 

20 In its simplest forms, ehe invention assumes 

certainty. That is, the value of each parameter, both of 
the systems and of the decision-making model that is 
needed to determine that solution, is kno%m. Boundary 
limits may be generated outside the model. 

25 From one aspect, the method of the invention 

conqprises the steps of finding a set of satisf iable 
boundary limits and developing in computer memory a 
multitude of feasible candidate strategies that satisfy 
the boundary limits, finding the values of outcomes of 

30 the risk- related activities for ^ach feasible candidate 

strategy under relevant conditions by storing and solving 
a model of the physical systems in computer memory, 
identifying candidate strategies from the outcomes of the 
risk- related activities and recording such outcomes in 

35 computer memory, and jointly applying in conqputer memory 
multiple optimization criteria to the outcc^s of the 
candidate strategies to aid decision-makers In selecting 
an implementable strategy. 
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From another aspect, the method of the invention 
caa^riseB the steps of finding a set of satisf iahle 
boundary limits and developing in conputer memory a 
multitude of feasible candidate Btrat:egie8 that satisfy 
5 the boundary limits by computer methods that use risk* 
limiting constraints, finding the values of outcomes of 
the risk- related activities for each feasible candidate 
strategy under relevant conditions by storing and solving 
a model of the physical systems in computer memory, 

10 identifying the candidate strategies from the out<:omea of 
the risk- related activities and recording such outcomes 
in computer memory, jointly applying multiple 
optimization criteria in computer memory to outccnives of 
all candidate strategies to exclude candidate strategies 

15 which are inferior to others under the criteria, and 
jointly applying multiple optimization <:riteria in 
computer memory to outcomes of remaining candidate 
strategies to aid decision -makers in selecting an 
implementable strategy. 

20 Prom yet another aspect, the method of the invention 

comprises the steps of storing the initial decision- 
making model in con5>uter memory (the model including the 
known values for each needed parameter of the physical 
systems) , producing a solution of the initial model in 

25 computer memory, deriving from the solution of the 

initial model the levels of outcomes for risk- related 
activities, and coxnparing these levels with boundary 
limits. 

If some preliminary boundary limits are violated, 
30 the next step of the invention is to modify the decision- 
making model so that it explicitly includes risk- limiting 
constraints that eliminate or minimize these violations. 
If there exists no "feasible" solution, herein defined as 
one in which all preliminary boxmdary limits are 
35 simultaneously met, at least some of the boundary limits 
are changed. After such a x:hange, the e3q>emded initial 
model with added risk- limiting constraints - based on 
both earlier solutions of the model and changed boundary 



wo 98/13776 



26 



PCr/US97/16446 



limits - ie solved again, and its outcoznee in risk- 
related activities are compared anew with the boundary 
limits. This process continues iteratively until a 
feasible solution is obtained. 
5 After that, a similar procedure is used ^o obtain a 

"satisfactory" solution, which meets both predetermined 
and discretionary constraints. 

In distinction from the customary mathematical 
programming. lAere initial model conatrainte are 

10 considered "sacred" and have to be fully met, this 

invention admits the possibility of inaccurate data being 
used in formation of almost any constraints. (The 
■'predetermined" constraints based on valid srules and 
regulations may be an exception.) After all, even if a 

15 model cannot find a feasible solution, real life as a 
rule still finds it, showing that the model is wrong. 
Therefore this invention does not necessarily battle 
infeasibilities, as done in standard mathematical 
programming or Robust Optimization, but instead prefers 

20 to analyze and change the constraints, whenever 

appropriate and necessary, to find r.ealistic feasible 
solutions. 

A key feature of the invention is that it aids in 
risk management under conditions of uncertainty. In such 

25 a situation, at least one parameter of a physical system 
has a plurality of different possible values. Prom this 
aspect of the invention, the method aids in risk 
management under uncertainty through the use of 
scenarios. A "scenario" is here defined as a 

30 respectively different combination of values of all 
needed systems parameters, with each combination 
including a variation in the value of some uncertain 
{i.e., multiple value) parameter. The decision model 
tinder uncertainty from the very beginning includes risk- 

35 limiting constraints that stem from predetermined 
boundary limits on risk- related activities. 

From this aspect of the invention, the method 
conprises the steps of storing a decision model of the 
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systems in computer memory , the model including ^he 
values for parameters of each physical system and at 
least one of the parameters having a plurality of 
possible values because of uncertainty; producing 
5 solutions of the model in -computer memory for a plurality 
of scenarios; and deriving trotn the solutions of the 
model the outcome levels for the rislc- related activities 
of each of the used scenarios. 

If, for some scenarios, there exists no feasible 

10 solution in which all boundary limits on risk-related 

activities are met simultaneously, the set of limits may 
be changed. After such a •change is made, the model is 
applied again with modified risk-limiting const^ainta 
that are based on changed limits. 

15 Prom yet another inportant aspect of the invention, 

the coniputer method aids in risk management under 
uncertainty by using scenarios via either single -scenario 
models or multiscenario models that have scenario 
submodels. (A single- scenario model may be considered as 

20 a special case of a multiscenario model, where it also is 
a single submodel.) Either way, the method of this 
invention is to develop candidate strategies - seta of 
niraierical values for outcomes of a relatively small 
number of key (i.e., strategic) variables that correspond 

25 to either the most inqportant or immediate and irrevocable 
decisions. In both instances, the method also focuses on 
an '•outcome matrix" - a three-dimensional array of 
outcomes of risk-related activities, systematically 
arranged in a "scenario versus strategy versus risk- 

30 related activity" format. (The number of dimensions may 
be greater than three.) A "regret matrix- derived from 
the "outcome matrix* plays a no less inqportant role. 

The con5)ressed, highly aggregated results 
accumulated in the outcome and regret matrices serve two 

35 main purposes. First, they allow easy generation of a 

major category of risk-limiting constraints and therefore 
help to form and modify the candidate strategies. 
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Second, they permit the etrategles to be analyzed and 
compared. 

When single -scenario models are used, the method of 
the invention generates strategies from a multitude of 
5 solutions of single -scenario models. It produces 

solutions that are both • feasible" (defined as meeting 
all predetermined boundary limits) and "satisfactory" (as 
escplained before, defined as meeting the "discretionary" 
limits that are imposed during the computational process 

10 of strategy formation) . The method then constructs the 
candidate strategies by classifying the values of 
strategic variables in each scenario solution into 
specific segments. Finally, it adapts the strategies to 
conditions of all scenarios and places the resulting 

15 outcomes into the corresponding cells of the outcome 
matrix. 

When multiscenario models are used, the method 
stores an initial multiscenario model of the systems in 
computer memory. The initial model includes parameters 
20 for each physical system, with at least one parameter 
being uncertain in having a plurality of different 
possible values. The multiscenario model conqprises 
interconnected scenario submodels for each scenario, each 
submodel being filled in computer memory by the values of 
25 parameters from the corresponding scenario. The 

multiscenario model produces a compromise solution that 
talces into account all scenario sukxnodels, as well as 
probabilities of scenarios, and includes a solution for 
each of these submodels. The method derives the outcome 
30 levels of the ris)c- related activities for the respective 
scenarios from solutions of the submodels. 

If the decision-makers are not satisfied with the 
outcome levels for some rislc- related activities under 
some scenarios, or want to find and explore some 
35 tradeoffs, the method of the invention adds 

"discretionary" risk- limiting constraints -to the initial 
multiscenario model and produces solutions -to the 
expanded model with these risk -limiting constraints that 
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move the solution in the desired direction. This 
addition of discretionary risk-limiting constraints is 
equivalent to overriding the original scenario 
probabilities, which are the least reliable element of 
5 input information. 

Prom another aspect of the invention » the 
discretionary risk-limiting constraints can be added to 
the initial multiscenario model consecutively, in ah 
iterative manner. The invention -thus allovrs the out<:ome 

10 levels of risk- related activities to be brought within 
acceptable limits on a step-by-step basis. 

From yet another aspect' of the invention, if the 
numbers of either scenarios or candidate strategies are 
too large and require an excessive con^utational effort » 

IS each of these numbers can be reduced to an acceptable 
level - either by clustering "individual" scenarios or 
strategies into "typical" scenarios or strategies, 
respectively, or by using some other basis for grouping 
scenarios, such as listing functional or statistical 

20 relationships. In distinction from the usual 

applications of clustering, where errors arising from 
incorrect composition of clusters are difficult to find 
and eliminate, this invention, when dealing with doubtful 
cases of assigning individual scenarios or strategies to 

25 groups, can "uncluster" them fully or partially and 

consider/ whenever necessary, either the "individual" 
entities instead of "typical" ones or reclustered 
"typical" entities. 

From yet another aspect of the invention, the method 

30 may also deal with a large number of scenarios and with 
the difficulty of solving the resulting large-scale 
optimization models. Instead of finding the solution of 
one large-scale multiscenario problem that covers S 
scenarios, the method in turn solves many smaller 

35 multidimensional models of a size that can be easily 

handled - e.g., S/IO of -ten -scenario models. This method 
of the invention conibines the advantages of solving smsai 
models with the possibility of finding compromise 
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Strategies, well adjusted to combinations of diverse 
conditions. It also introduces new limiting constraints, 
which are based on the values of strategic variables in 
solutions of multiscenario problems formed for groups of 
5 worst scenarios • This method is a relatively complex 
aspect of the invention, but is also particularly 
effective for dealing with very large and canapl±<zatBd 
problems. 

By changing the conqposition of scenario groups 
10 covered by each multiscenario model and by fine-tuning 
the discretionary risk* limiting constraints in^osed on 
the models, the invention permits finding good candidate 
strategies in a mere fraction of the time needed for such 
confutations in existing state-of-the-art methods. 
15 As for methods that deal with selection of the 

candidate strategy to be implemented in physical systems, 
another aspect of the invention is to provide three new 
synthetic regret -based methods. 

Prom yet another aspect of the invention, these 
20 methods are presented in an analytical and graphic 
framework of newly developed ^efficient frontier 
(regret) . " 

Another aspect of the invention involves 
simultaneous use of multiple weights, or conversion 

25 coefficients, to sin^lify either the outcome and regret 
matrices or the process of decision-making. Sin^lifying 
the outcome or regret matrix may mean it is converted 
from three or more dimensions to two dimensions, or that 
qualitative parameters are converted into quaujtitative 

30 A^lues. Sinqplifying the process of decision- making may 

mean that multiple criteria are amalgamated into a single 
criterion, or that the opinions of different members of 
the decision- making group are made commensurable. 
Similarly to other mathematical methods, such 

35 simplification is often not reliable. Therefore, as in 
the stage of constructing the strategy, conversion 
techniques are primarily used to screen out bad and risky 
strategies, rather than to select "the best" strategy. 
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Moreover, the invention combines the conversion eature 
with imposing (at the strategy construction stage) rislc- 
limiting constraints, which stabilize and moderate the 
decisions. 

5 The invention may be more fully understood from the 

following detailed description of several specific 
exan9>le8, considered in the light of the accompanying 
drawings and the appended claims. In flow charts « 
parenthetical numbers are keyed to similarly numbered 

10 mathematical models in the detailed description text« 
Brief DeacriPtion of the Drawings 

PIG. 1 is a flow chart of a computer method, in 
accordance with the invention, for managing risk under 
conditions of certainty; 

15 PIG. 2 is part of a flow chart of a con^ter method, 

in accordance with the invention, for managing risk under 
conditions of \incertainty, using single -scenario 
models; 

FIG. 3 is a continuation of the flow charts in PIG - 
.20 2 and FIG. 5; 

FIG. 4 is a flow chart of a <:omputer method, in 
accordance with the invention, for msmaging risk under 
uncertainty, using a single multiscenario LP model; 
FIG. 5 is a flow chart of a computer method, in 
25 accordance with the invention, for mamaging risk under 

uncertainty, using single -scenario UP models jointly wi-th 
clustering individual scenarios into typical scenarios; 

PIG. 6 is a flow chart of a computer method, in 
accordance with the invention, for managing risk under 
30 uncertainty with the use of multiple multiscenario LP 
models; 

PIG. 6A an extension of the flow chart in FIG. 6; 
FIG. 6B is another extension of the flow chart in 
FIG. 6; 

35 FIG. 7 illustrates a two-dimensional payoff matrix 

for a four- strategy, six scenario problem; 

FIG. 8 illustrates a two-dimensional regret matrix 
for the same four -strategy, six scenario problem; 
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FX6. 9 Is a decision graph (OPR efficient frontier) 
for tne OPR method of <:azi^arlng the candidate strategies; 

PIG- 10 is a decision graph (PIR efficient frontier) 
for the PIR method of comparing the candidate strategies; 

FIG. 11 is a decision graph (MIRR efficient 
frontier) for ^he MIRR method of comparing the candidate 
strategies; and 

FIG. 12 is an extension of the flow charts of FIGS. 
3, 4, and €B illustrating selection of strategies. 

It is believed useful, at this point, to provide an 
overall view of the invention before discussing details 
and the figures. 

This invention is a hybrid, an ensenible of two major 
directions in methods of decision-ma)cing, <hose of OR/MS 
and those of DS. The two technologies con5>lement each 
other and are critical to each other's Buccees. The 
invention uses the inproved tools of DS for the selection 
of the strategy, but it also uses OR/MS optimisation 
models for formation of candidate strategies and -for 
filling in the outcome matrix. Moreover, the 
optimization models are used in such a way that, instead 
if increasing ris)c, they reduce it. Realizing that 
mathematical models are never completely reliable in real 
world decision -ma)cing, and that optimization models are 
especially vulnerable because of their introduction of an 
additional risk component of their own, the invention 
uses optimization primarily to screen out bad and ris)cy 
strategies, rather than to construct and select "the 
best" strategy. This single but powerful technique helps 
to deal successfully with very coii^li*cat:ed, previously 
unsolvable problems. But this approach also turns out to 
be the essence, indeed, the very goal of risk management. 

This invention uses techniques of prior art in OR/MS 
and DS, but uses them in different framework and €or 
different purposes • 

First, the primary tool of this invention is the 
explicit addition to the main model (or a formula or a 



wo 98^3776 



PCTA)i597/16446 



method) of multiple risk- limiting constraints. Mosc of 
the explicit additional conetrainte are of a "cutting 
plane" type, well-known in the prior art. These 
constraints were first introduced into mathematical 
5 programming by Ralph E. Gomory in "Outline of an 

algorithm for integer solutions to linear progrrams," 
Bulletin of the American Mathematical Society, v. 64 
(195B), pp- 275-278. Up to the present, cutting 
hyperplane constraints have been used to cut off non- 

10 integer portions of the region of feasible solutions and 
some other similar purposes, but in this Invention 
similar constraints Instead cut off portions of the 
feasible solution region that are too rislcy, where the 
unacceptable outcomes are located, leaving for further 

15 optimization only a trxmcated portion of the legion where 
the risks are acceptable. The Invention may also use 
such forms of risk- limiting constraints as fines and 
penalties that perform the same role, namely, they limit 
the possibility of undesirable outcomes, but in a weaker 

20 (non- mandatory) form. More specifically, while the risk- 
limiting constraints in the form of equations and 
inequalities set the absolute limits (boundaries) of 
outcomes in risk- related activities that cannot be 
exceeded in the LP model solution, the constraints in the 

25 form of fines and penalties adversely affect the movement 
toward or over the boundaries, but do not categorically 
forbid it. 

The risk- limiting constraints, including fines and 
penalties, can be added to the LP model in any order and 

30 way desired. Xf they are added Iteratively, one by one 
or in groups, starting from different out<:ome levels for 
various risk- related activities, the resulting truncated 
feasible regions may differ, as differ the solutions of 
the evolving models. This procedure leads to the 

35 generation of different candidate strategies, with 
diverse acceptable tradeoff, conibinations of outcome 
levels in various risk-related activities. 
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Optimization models are full of constrainte of both 
types, but these predominantly are either constraints on 
the available capacities and other resources, or 
constraints on meeting the given production targets, or 
5 "balancing** equations that balance inputs and outputs in 
space and in time. In contrast, risk* limiting 
constraints are constraints on results. 

In turn, risk- limiting constraints on results can be 
classified in two categories; both can be of 4:he 

10 equation- Inequality type and of the fine-penalty type. 

The first category is "predetermined" constraints that, 
for instance, prohibit too great amount of waste water to 
be generated at a plant. These constraints are based on 
established norms, laws, and regulations that are known 

15 in advance. The constraints of the second <:ategory are 
not known in advance: they are based on <:onqputatlonal 
results attained for the models. 

In spite of the obvious advantages of adding 
es^iicit risk- limiting constraints of the second 

20 category, only one source is known to apply them, 

although in their simplest embryonic form. In "The uses 
of previous experience in reaching statistical 
decisions," Annals of Mathematical Statistics, 23, 396- 
407, 1952, Hodges and Lehmann suggest finding a decision 

25 that minimizes the maximum possible risk; let this 

maximum risk value be C. Then the level of "maximum 
tolerable risk" is chosen that is lower than C. The use 
of maximum tolerable risk as an auxiliary criterion in 
selecting a best strategy is tantamount ^o imposing a 

30 single and slxqple risk- limiting constraint of the second 
category . 

This Invention not only widely usee multiple risk- 
limiting constraints of the second category, it also 
introduces a new type of such constraints, namely, 
35 constraints on results of individual risk- related 
activities under individual scenarios or groups of 
scenarios. The constraints are mostly discretionary and 
allow the user to find and analyze best tradeoffs. The 
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constraints o£ this category can be easily applied only 
in multiscenario optimization models specific for this 
Invention that come next. This explains why such 
constraints have not been used previously. 
5 Second, this invention uses multistage, 

multiscenario LP models that include the risk- limiting 
constraints above. Similar but simpler models are 
utilized in stochastic programming. However, the 
multiscenario models make a novel and extremely effective 

10 combination with risk-limiting constraints. 

Third* this invention constructs alternative 
strategies, tests ^heir performance under different 
scenarios, evaluates the risk associated with each 
strategy, and tailors each candidate strategy to the risk 

15 attitude and personal preferences of DMs. 

For that purposes, this invention introduces two 
concepts, "strategic variables" and "candidate strategy." 
"Strategic variables" originate from optimization models, 
such as LP models described above; they usually are a 

20 small subset of variables selected from the whole set of 
the variables of the model. Strategic variables are ^he 
key decision variables, especially 4:he variables 
associated with immediate and irrevocable decisi<ms, that 
should keep their values stable or even constant tinder 

25 different conditions (scenarios) . For instance, in ^che 
distribution problem described above, the subset of 
strategic variables may be limited ^o the production 
targets at all industrial enterprises of the system. All 
other quantities to be determined, such as inventory 

30 levels or transportation flows, correspond to less 

ixnportant, non- strategic variables of the model. If <he 
number of strategic variables in a problem is K, then a 
set of values of all K strategic variables contained 
within the same segment defines a "candidate strategy." 

35 The strategic variables used in the present 

invention are partly similar to the "first -stage 
variables" in two- stage and multietage stochastic LP 
models; see, for instance « Chapter 16 in Harvey M. 
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Wagner, ■ Principles of operations Research," Prentice- 
Hall, Inc.. Englewood Cliffs, NJ, 19€9, pp. 657-€€8 and 
678-684. However, strategic variables are different from 
the first stage variables. They do not necessarily 
5 belong to the first stage. Furthermore, there is no 
concept in OR/MS similar to the "candidate strategy," 
where the strategy is defined not by the total solution 
of an LP model, but rather by the values of a small 
subset of a model's variables. 
.0 Fourth, this invention focuses on a relatively small 

niunber of strategic variables/ instead of all variables 
of the computational model used, and scenario submodels 
of the multiscenario model are connected only through 
values of strategic variables. Combined with an 
L5 understanding that a good solution of a multiscenario 

model must also be good under the worst scenarios, this 
feature brings about unexpected but crucial computational 
benefits. This invention offers a novel, highly 
efficient and quickly converging deccHoposition method, in 
20 which the boiindary limits on the values of strategic 

variables are derived from solutions of worst scenario 
s\ibmodel8, and then extended to all other scenarios. 

Fifth, this invention introduces novel, synthetic 
regret -based methods, or optimality criteria of DS, for 
25 comparing candidate strategies. They are accompanied 

also by a novel, "regret-based efficient frontier," which 
is developed both in analytical and graphic foann. 

The definition of efficient frontier in Portfolio 
Theory has been given earlier; in this invention, 
30 efficient frontier is defined differently, as the set of 
strategies that are the best for the €\xll range of 
eaqiectations, from the optimistic to the pessimistic end. 

It should be pointed out that these methods can, as 
a practical matter, be executed only on a conputer, 
35 because they are used within a complex computer-based 
framework with multiple interconnected criteria and 
decision graphs, and in some (especially financial) 
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applications may have to be repeated a large niunber of 
times. 

Sixth, this invention combines the multitude of 
jointly used diverse DS criteria by means of conversion 
5 coefficients, or weights. Obviously, the values of these 
<:oeff icients can be known at best approximately. 
Therefore the joint use of many criteria is considered 
not as a final and conclusive step, but just as one of 
many ways to analyze and compare strategies. 

10 Weighting diverse optimality criteria in 

multlcriterion models to merge them into a single 
criterion model, such as LP, is a technique well known in 
state-of-the-art methods. Weighting such disparate model 
results as, e.g., %»8te water and financial loss, may 

15 also have been used before. This invention introduces a 
new element to the weighting approach, whereby several 
values of the weights are used simultaneously, which is 
achieved by increasing the number of scenarios. 

This invention has several major conponents. It: 

20 (a) Adds risk- limiting constraints (which include 

equations, inequalities, fines, and penalties) -to the 
main decision-making formula or model for decision-making 
both under certainty and uncertainty; 

(b) Offers a method of risk protection and 

25 management by imposing risk-limiting constraints that are 
not only more reliable but also simpler than those 
generated by state-of-the-art methods and therefore can 
be added to the main model or formula by people not 
skillful in the art of modeling; 

30 (c) Develops cwpromise ^candidate strategies on the 

basis of either multiple single- scenario, or single 
multiscenario, or multiple multi scenario mathematical 
programming models, all with additional risk- limiting 
constraints; 

35 (d) Solves extremely large multiscenario 

mathematical programming models by a novel decomposition 
method that efficiently uses the information derived trcm 
solutions of submodels for worst scenarios; 
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(e) Develops multiple scenarios and multiple 
candidate strategies; groups and clusters scenarios and 
strategies, whenever this is necessary because of the 
coav>utational effort involved; and declusters, reclueters 

5 and regroups the scenarios and strategies whenever these 
operations can improve the analytic and decision support 
capabilities of the method; 

(f ) Presents the "strategy versus scenario" results 
in the format of two multidimensional matrices, an 

10 outcome matrix and a regret matrix, and then applies 

different methods and criteria of DS, including special 
methods of this invention, to help, first, in screening 
out the worst strategies, and second, in selecting a 
strategy; 

15 (g) Uses special synthetic methods for selecting a 

strategy based upon coxi?>aring the potential reduction of 
risk (regret) with the sacrifices needed to achieve that 
risk reduction, as well as upon weighting the outcomes 
achieved under different criteria; 

20 (h) Compares the candidate strategies in an 

analytical and graphical framework of "efficient frontier 
(regret) , " which not only provides excellent help in 
decision- making, but also minimizes the negative intact 
of uncertainty; 

25 (i) Simultaneously applies several values of weights 

to combine: diverse optimality criteria; model outcomes 
in different risk-related activities; opinions and 
attitudes of different decision -makers; and different 
qualitative characteristics . 

30 In FZ6. It a flow chart illustrating application of 

the invention to the management of risk under conditions 
of certainty begins with step 11, in which an initial 
symbolic model of the physical systems of interest is 
stored in conqputer memory. From there, the method moves 

35 through step 13, in which parameter values are entered in 
memory, to step 15, in which boundary limits which are 
not to be exceeded are entered in memory. Prom step 15, 
the method moves through step 17, in which a 
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computational model of the physical systems is generated 
in computer memory « to step 19, in which the 
coinputatlonal model Is solved. Prom step 19^ the method 
moves through step 21, where I is the number of t:imes the 
5 computational model has been solved and is initially set 
equal to 1, to step 23, where outcome levels of the 
various activities of the physical systems are obtained. 

In PIG. 1, the illustrated method moves from step 23 
to decision point 25, where the outCOTie levels from step 

10 23 are cosqpared with the boundary limits entered in step 
15. If the answer at decision point 25 is yes (i.e., all 
outcome levels are within boundary limi-ts) , the method 
moves to and terminates in step 27, in which the boundary 
limits are implemented in the physical system or systems. 

15 xf the amswer at decision point 25 is no (i.e., at least 
one outcome level is outside of the boundary limits), the 
method moves to decision point 29, where it is determined 
whether or not I is still equal to 1. If the ansvrer at 
decision point 29 is yes (i.e., the model has only been 

20 solved once), the method moves to step 31, in which risk- 
limiting constraints triggered by the solution in step 19 
are added to the cooqputational model. The next step is 
then step 33, in which the esqpanded model is solved. 
Following step 33, the method moves to step 35, in which 

25 I is incremented by 1 to indicate that the model has been 
solved another time (e.g., for a second solution, 1-2). 
The method then moves baclc to step 23, where outcome 
levels are obtained. 

After step 23, the method moves again to decision 

30 point 25, where it terminates in step 27 if the answer is 
yes and moves to decision point 29 if the answer is no. 
If the answer at decision point 29 is no (as it will be 
if I is greater than 1) , there will be no feasible 
solution unless the boundary limits are modified, so the 

35 method moves to step 37, in which one or more of the 
boundary limits are modified. Following step 37, the 
method moves to step 39, in which the risk- limiting 
constraints added to the model at step 31 are modified. 
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Prom step 39 , the method proceeds to step 33 « in which 
the expanded model is solved again. 

As explained in the "Background" section o£ this 
specification, decision- making problems smd situations 
5 usually are very complicated. Unable to deal properly 
with all complexities, humans make decisions following 
some approximate "^rules of thumb" that may prodtice both 
good and bad decisions. The most glaring coonmon flaw of 
most "rule of thumb" decisions is insufficient attention 

10 to risk. 

A main premise of this invention is that, in order 
to protect and insure physical systems from risk, 
decision-makers (DMs) must es^licitly limit the outcome 
levels of risk-related activities, usually in the form of 

15 risk-limiting constraints. This brings about two 
important consequences. First, the addition of 
constraints transforms even single models into more 
complicated versions, which caxuot be solved in 
rea8onsUt>le time without cocnputers. Second, the state- of- 

20 the -art optimization models currently used in these 
complicated decision-making situations, when applied 
without due attention to risk, are insufficient and often 
counterproductive, because any optimization model adds a 
substantial risk con^onent of its own. This invention 

25 provides the techniques necessary for overcoming such 
difficulties. 

As a rule, both these consequences are true ev^n in 
the relatively simple case of decision-making under 
certainty. Consider a physical system under certainty, 

30 where DMs know the exact single value of each system 

parameter that is needed to determine the future behavior 
of the system under different conditions and, therefore, 
to make a decision. The system performs interrelated 
activities. Some of them are risk-related, that is, the 

35 outcome levels for these activities may fall outside some 
predetermined boundary limits, and that causes -concern 
for DMs. The decisions to be made about the system 
specify the planned levels for some, but not all, 
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activities, such as the utilization of limited or ^costly 
resources, or the outputs o£ some products. Because of 
the complexity of both internal and external connections, 
relationships, and reactions of the system, the outcomes 
5 of other activities. Including the risk- related ones, are 
not obvious. These other outcomes are determined by a 
special model of the system. The unique set of outcome 
values for all activities of the system, defined by the 
model on the basis of the decisions to be made, is called 

10 "a solution. The model also includes a set of criteria 
that determine which solution is "the best," that is, 
which one to choose for implementation among many 
candidate decisions. The model thus not only describes 
the physical system, but also assists in decision-making. 

15 Suppose that, using the model of the system, DMs 

made a decision and determined all outcomes coming from 
that decision. Suppose also that, after coxnparing the 
outcome levels for risk- related activities with 'the 
predetermined boundary limits, DMs discovered that some 

20 outcomes fall outside these limits. As a rule, as 

indicated in the "Background" section, the only way to 
guarantee meeting all limits is then to explicitly append 
to the initial model the additional risk- limiting 
constraints, usually in the form of ^equations and 

2 5 Inequal 1 1 ies . 

This rule has an exception. Suppose ^hat we have a 
simple system, such as choosing an appliance like a 
washing machine. The model of the system is also 
extremely sinqple: it is just a list of "attribute 

30 packages" that assigns -to each type and make of machine 
the relevant values of performance attributes, such as 
features, quality, price, and service levels. The 
selection procedure is just coxnpiring a list of machines 
that have the desired attributes (^<:apacity no less than 

35 , quality no less than etc.") and then choosing 

from that list the lowest price machine. Let <he risk- 
related activity be the monthly cost of energy, which was 
not. included among the attributes, and the DMs want this 
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cost not to exceed x dollars. DMs can avoid here ueing 
the model with an explicitly added risk- limiting 
constraint in the form of an equation or inequality, such 
as "monthly energy cost no more than $x, " in two ways: 
5 either this item is added to the attributes and then -the 
list of machines is further purged by removing all 
machines that do not meet this criterion; or, if a 
machine's monthly energy bill is $y and y > x, -the 
selection is performed not on the pri<:e of the machine, 

10 but on the total of price plus fine, where the fine may 
depend on the value of (y • x) . 

If the model of the system is more <;omplex than just 
a list of "attribute packages," so that the risk-limiting 
constraints in the form of equations or inequalities must 

15 be added to the model explicitly, the resulting model 

becomes one of a "constrained optimization" mathematical 
programming type, which as a rule can be solved in 
reasonable time only on a computer. There are various 
versions of mathematical programndng, such as linear 

20 programming (LP), non- linear programming, integer 

programming and so on. For sixtiplicity, this discussion 
will not go beyond LP, but the arguments and conclusions 
dp not depend on which type of optimization model is 
used. 

25 Suppose that, prior to addition of risk- limiting 

equations and inequalities, the decision- support LP model 
is described, in matrix notation, by (1), whioh contains 
expressions (1.0) through (1.2): 

Find a (n x 1) vector x to 
30 Maximize ZD « C*^X il.O) 

subject to AX - B (1.1) 
and X 2 0. (1.2) 

Superscript represents a matrix transposition 
operation. Model (1) comprises: a vector X of variables 
35 x(q) (q - 1, . n) that define the levels of 

activities q of the physical system; ro equations (1.1) 
iiqposed on limit d resources, product outputs, and 
interconnections between activities q; n non- negativity 
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conditions (1.2) on variables x[q] ; and an objective 
function ZD - a scalar product of tvo vectors C and X of 
dimension n defined by (l.O) • The value of ZD is the 
difference between benefits and costs; it is to be 
5 maximized. 

All coo^ponents of vector C » (c(q]) and (m x l) 
vector B » (b[p]) and of (m x n) matrix A » [a(p] (q] ] (p 
-1, ... , m; q » 1, . . . , n) are calculated on the basis 
of values of parameters of the physical system and stored 

10 in the computer memory. Variables x(q] correspond to 

levels of activities q (^hat is, scale of operations or 
output of products or use of resources) of n 
interconnected activities q. Equations (1.1) are m 
linear equations with coefficients a(p] [q] on n variables 

15 x[q] . These equations have right-hand sides b[p] . 

Equations <1.1) connect levels of operations or specify 
the output of products or limit the use of resources -for 
activities q. In the exceptionally simple case of an 
appliance described above, there are no equations (1.1), 

20 so m « 0. Objective function (1.0) is a linear function 
of varieO^les x(q] • 

As mentioned before, the last R of n activities of 
the physical system (j « n - R 4- 1, ... , n) are risk- 
related activities that are of special conceam to the 

25 DMs. The values of the corresponding variables x{jl 

determine levels of these activities, such as the outputs 
of some dangerous by-products. For simplicity, here and 
in subsequent models (3) and (4) assiime that the 
predetermined boxuidary limits for these R activities are 

30 one-sided, that is, the outcome for an activity is 

limited by only one upper boundary - say, it should be no 
greater than this boundary limit. 

Again, if DMs have solved model (1) and have 
obtained the values of outcomes for all risk-related 

35 activities, and if all outcomes are no greater than the 

corresponding boxmdary limits, then the decision baeed on 
that solution is acceptable to DMs, and the con^tations 
can be stopped. However, if at least one outcome exceeds 
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the boundary limit £or the corresponding activities « thien 
the model (1> must be expanded into the model (Xh) , whi<di 
contains esqpressions (l.OA) through (1.5A) : 

Find (n X 1) vector X and (R x l) vectors W auid Y 

5 to 

Maximize ZD = c'^X - F^Y (l.OA) 
subject to AX a B, (I.IA) 
Xj^ 1- W - y BL. (1.2A) 
X 0, (1.3A) 

10 W j: o/ (1-4A) 

and Y * 0, U.5A) 

where: Xj^ is an (R x 1) vector of outcomes (where Xj^ is a 
subset of the set x) for risk- related activities, ndiich 
are the last R of n activities; BL is an <R x 1) vector 
15 of boundary limits for the same set of activities; w and 
y are, respectively, (R x l) vectors of zero priced 
"slack variables" and highly fined "artificial variables" 
that transform risk- limiting inequalities into equations; 
and P is a (1 x R) vector of fines f Iql (q - N - R ♦ 1, 
20 ... , n) that is imposed on the excessive portions y[ql 
of outcomes that exceed the boundary limits BL*{q] <q - n 
- R ♦ 1, ... , n) . Equations (1.2A) are R risk- limiting 
inequalities transformed into equalities by adding slack 
variables W and subtracting artificial variables Y, while 
25 (1.4A) and (1.5A) are non- negativity conditions for ^hese 
variables. The transformation of model (1) into model 
(lA) can be performed either all at once for all R risk- 
related activities, or iteratively, adding .equations 
(1.2A) only for those activities where the outcomes in 
30 the previous solution exceed the boundary limits. 

It is important to notice that, even if model (1) is 
single enough to be solved manually (for instance, if the 
number m of "structural" equations (1.1} that describe 
the physical system equals zero), model (lA) still 
35 becomes an LP model, suf f iciently <:omplicated to be 
solved in a reasonable time only on a computer. 

It is also iii«>ortant that ms cannot be sure that 
the values of the predetermined boundary limits BL[ql (q 
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oD-R^^l, .*« #n) are correct or correspond to real 
capabilities o£ the physical system, especially if these 
values do not allow for obtaining a feasible solution of 
model (lA) , that ie, a solution where all outcomes are 
5 within the boimdary limits and therefore Y « 0« 

Therefore DNs' actions may include modifying the values 
of these limits. 

The method of DMs* behavior in finding the solution 
for a physical system under <:ertainty *can easily be 
10 followed on the flow chart of FIG. l. Af^er formulating 
the initial symbolic model (1) and storing it in -coinputer 
memory together with the values of the relevant 
parameters and boundsury limits, tmbz 

generate the computational model (1) - ste{>8 11 - 

15 17; 

solve model (1) by any appropriate means and, if 
model (1) is an LP model, by any appropriate computer LP 
system - step 19; 

record in the computer memory ^he values of outcomes 
20 obtained from the optimal solution of model (1) for jeach 
of R risk-related activities - step 23; 

compare the outcome with the boundary limit for each 
risk- related activity and decide whether the obtained set 
of outcomes is fully acceptable -for all such activities - 
25 step 25; 

if this set of outcomes is fully acceptable, the 
method is com^jleted - step 27; 

otherwise/ DMs add to model (1) risk- limiting 
constraints, transform it into model (lA) - step 31, and 
30 solve model (lA) 1^ any appropriate computer LP system - 
step 33; 

again, conpare the outcome with the boundary limit 
for each risk -related activity and decide whether the 
obtained set of outcomes is fully acceptable for all such 
35 activities - step 25; 

if the set of outcomes is acceptable, the method is 
crapleted - step 27; 
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otherwise « modify the boundary limits and the 
corresponding risk- limiting constraints and solve the 
resulting model (lA) - steps 37, 39, and 33; and 

continue the procedure of modification of boundary 
5 limits until an acceptable solution of model (lA) is 

foimd, which completes the method - steps 37, 39, and 33. 

In PIG- 2, a flow chart illustrates a method 
applying the invention to the management of risk under 
conditions of uncertainty by using SSLP (single -scenario 
10 linear programming) models* The method begins with step 
41, in which an initial symbolic SSLP model (3) is stored 
in computer memory. Next, the method moves -to step 43, 
in which the parameter values *or S different scenarios 
are entered into coaputer memory. After step 43, the 
15 method moves to step 45, in which any predetermined 
boundary limits for the various outcome values are 
entered. Following that, step 47 enters value segments 
for K different strategic variables. Step 49 forms S 
scenarios, each made up of a different -coiribination of the 
20 certain parameters with a value for each uncertain 

parameter. After step 49 comes step 51, whioh sets H 
equal to 1 to indicate the index of the first of the S 
different scenarios, where H is the index number of the 
scenario being tested. 
25 in PIG. 2, the illustrated method next moves <o step 

53, where a conputational model (3) is generated in 
con?>uter memory for scenario H, After that oomes step 
55, where the computational model for scenario H is 
solved, and step 57, which obtains outcome levels from 
30 the solution of scenario H.. Following step 57, a 

decision point 59 determines whether or not the solution 
is feasible (i.e., whether all its outcome levels are 
simultaneously achievable). If the answer is yes, the 
method moves to decision point €1, which determines 
35 whether or not the outcomea are all satisfactory. If the 
answer is no, the method moves to step «3, which adds or 
modifies discretionary constraints to the model in 
computer memory to limit unsatisfactory outcomes. Next, 
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Step 65 sets H equal to 1 to indicate that the method 
begins again with scenario 1. The method returns then to 
step 53 to generate a modified or expanded cosputational 
model for scenario H « 1. 
5 If the answer at decision point 59 is no, the method 

moves to step 67 to modify the predetermined boundary 
limits. Next, the method moves to step ^9, modifying any 
previously entered risk- limiting constraints that are 
included into model (3) in computer memory. 

10 If the answer at decision point 61 is yes« the 

method moves to decision point 71 to determine whether H 
is equal to S, where S is the total number of acenarios 
being tested. If H is not equal to S (that is, less than 
S) , the method moves to step 73, where H is incremented 

15 by 1, and back to step 53 to generate the <:oinputational 
model (3). 

If the answer at decision point 71 is y^s, the 
method moves to the flow chart in FIG. 3, whi^h begins 
with step 74 and moves to step 75 to find ^he valiies of K 

20 strategic variables in S scenario solutions, where K is 
the total number of strategic variables. In FIG. 3, the 
method then moves to step 77, which classifies values of 
K strategic variables in S solutions into predetermined 
segments. Prom step 77, the method moves first to etep 

25 79, which defines NST candidate strategies, where NST is 
the total number of candidate strategies, and then to 
step 81, which sets I and J both .equal to zero, where I 
is the index of the strategy being tested (the range of I 
is from 1 to NST) and J is a count of ^he strategies 

30 which have already been tested and accepted. From there, 
the method moves to step S3, which sets the acenario 
index H equal to zero. 

From step 83 in FIG. 3, the method moves first to 
step 85, which increments I by 1, and then to step 87, 

35 which formulates strategy definition constraints for 
strategy I. The next step in the method is step 89, 
which increments H by 1. Prom step 89, the method moves 
first to step 91, which adds a model (3) for scenario H 
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to form a computational model (3A) , axid then to step 93, 
which solves this inodel for scenario H. After step 93 
comes decision point 95, where it is determined whether 
or not a solution to the computational model (3A) is 
5 feasible in the sense that it meets all predetermined 
hoiindary limits on risk-xelated activities. 

If the answer at decision point 95 is no, the method 
moves to decision point 97, where it is determined 
%rtiether H is equal to S- If the answer is no, the method 
10 returns to step 83. If the answer is yes, the method 
proceeds to decision point 99, where it is determined 
whether I is equal to NST. If the answer is no, the 
method again returns to step 83. If the answer is yes, 
the method proceeds first to step 101, whi^ is the 
15 strategy, selection process illustrated in PIG . 12 , and 
then to final step 103, which implements the selected 
strategy in the physical system or systems. 

If the ans%^r at decision point 95 in FIG, 3 is yes, 
the method moves to decision point 105, which determines 
20 *rtiether H is -equal to S. If the answer at decision point 
105 is no, the method returns to step 89 . If the answer 
is yes, the method moves to step 107, vrtiich increments J 
by 1. From step 107, the method moves first to step 109, 
which defines candidate strategy J, and then to step 111, 
25 which fills "strategy J" columns of a matrix of outcomes. 
From step ill, the method proceeds to decision point 113, 
which determines whether I is equal to NST. If the 
answer is no, the method returns to step 83. If the 
answer is yes, the method goes on to step 101 and, ^rom 
30 there, to final step 103. 

Two difficulties (i.e., the complication of a 
decision- support model and the insufficiency of <:urrently 
used approaches to applying optimization models) involved 
in making decisions about physical systems are outlined 
35 in the beginning of this detailed description. They are 
even more important if the system is under uncertainty, 
that is, when some parameters of -the system are not known 
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exactly and therefore may be represented by not one value 
but rather by a multitude of different possible values. 

Suppose we are to make decisions about a physical 
system with multiple parameters that performs n 
5 interrelated activities* Some of these activities are 
use of resources, manufacturing of products, storage or 
transportation activities, and so on. Let R <n > R k l) 
of these n activities be potentially dangerous, risk- 
related activities. The decisions to be made are about 

10 the planned levels of these n activities; each of these 

levels must be non-negative. The decisions will be based 
on the values of multiple parameters of the physical 
system. X>ue 'to uncertainty, G (<3 & 1) of these 
parameters have g[qj > 1 values <q - G) ; .each of 

15 the . remaining parameters has only one value. The maximum 
number MS of all respectively different combinations of 
possible values of all parameters of the physi<:al system 
is given by (2) : 

6 

20 MS « IT g[il ; <2] 

i-1 

where the symbol ir is used to denote "product of . " 

Bach of these respectively different combinations of 
possible values of all parameters of the system we will 

25 call a "scenario." Then we can f orm *S (l < S s MS) 

scenarios and store these scenarios in conputer memory. 

To construct candidate strategies, the method of 
this invention uses an LP model, which in the general 
case is a multiscenario linear programming model (MSLP) 

30 that has LS scenarios, where LS may vary from one to S. 
If LS < S, then, to cover all S scenarios of the 
problem, the method runs multiple MSLP models, each 
having up to LS scenario submodels. The greater the 
portion of S covered by LS, the more effective is the 

35 method. However, sometimes, due to the size of the LP 
problem, LS has to be small. 

In an extreme case, when IiS « 1, MSLP becomes a 
single- scenario LP (SSLP - single- scenario linear 
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programming) model. Suppose that, given one scenario H« 
the conqputer generates for that scenario H an SSLP model. 
This model has, however, a major distinction £rom model 
(1) . As mentioned before, the last R of n activities q 
5 of the physical system (q «» n - R + 1, — , n) are risk- 
related activities that are of special concern to DMs; 
the values of the corresponding variables t j J [H] define 
levels of these activities « that is, for instaui^e, the 
output of some dangerous by-products under scenario H. 
10 Decision-malcing under uncertainty involves more 

operations than decision -making under <:ertainty, and it 
is very likely that, at some operation, rlsk-ldLmiting 
constraints for predetermined boundary limits in R risk- 
related activities would be included in the decision 
15 model. Model (3), which, from the very beginning > 

includes these constraints and thus is similar to model 
(lA) , should therefore serv^e as the principal model, both 
symbolic and computational. 

For each scenario H, H » 1,...S, model (3) in matrix 
20 notation contains expressions (3.0) through (3.5): 
Find (n x l) vector XtHl 
and (R X 1) vectors WtHl and Y[H1 to 
^faxlmize ZM 1.0 x ZSIH] 
- 1.0 X CtHl'^ XIHl - F*^ YCHl (3.0) 
25 subject to A[H1X[H] » Btfl], (3.1) 

Xr[H] ♦ W(H1 - Ym - BL, (3.2) 
X(H1 2 0, (3.3) 
W[H} It 0, (3.4) 
and Y(W] K 0. (3.5) 

30 Superscript ^ represents a matrix transposition 

operation. Model (3) comprises a column vector X[H1 of 
variables X(q] [Hi (q « li..., n) that define the level of 
activities q of the system under scenario H; Xj^[K] is an 
(R X 1) vector of outcomes (where Xj^ is a subset of the 
35 set X) for risk- related activities, which are the last R 
of n activities? BL is an (R x 1) vector of boundary 
limits for the same set of activities, and these limits 
are invariant of scenarios; W[H) and Y(H) are. 
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respectively, (R x 1) vectors of 2ero priced "slack 
variables « and highly fined •artificial variables" that 
transform risk- limiting inequalities into equations; and 
P is an <1 X R) vector of fines F[ql (q « n - R ♦ 1, ... 
5 ^ n) , also invarismt of scenarios, that is imposed on the 
excessive portions y[q] [H] of outcomes that exceed the 
boundary limits BL[ql (q - n - R ♦ X, . . . , n) . An 
objective function ZStH] is defined by (3.0). The value 
of the objective function ZS[H], which is the difference 

10 between benefits and costs, is to be maximized » 
Equations (3.2) are R risk- limiting inequalities 
transformed into equations by adding slack variables W[H] 
and siibtracting artificial variables Y[H}, while (3.4) 
and (3 ,5) are non- negativity conditions for these 

15 variables, inequalities (3.3) guarantee non -negativity 
conditions for variables X[H) . 

All components of (1 x n) vector C[H] » (c[ql [H] ) 
and of (m X 1) vector B(H] - (b[p] [H] ) and o€ (m x n) 
matrix A IH J - [a[p][qlIHll (p - 1,.., m; q n) are 

20 calculated amd stored in the computer memory on the basis 
of values of all parameters of the physical system in 
scenario H. Variables xCq] [H] correspond to decisions 
that, under scenario H, define levels of activities (that 
is, scale of operations or output of products or use of 

25 resources) in n interconnected activities q. Equations 
(3.1) are m linear equations with coefficients a[p] [q] [H] 
on n variables xtql (Hi- These equations have right-hand 
sides b[pl{Hl. Equations (3.1) connect levels of 
operations or limit the output of products and use of 

30 resources in activities q. Objective function (3.0) is a 
linear function of variables x(p] [HJ • 

In addition to developing scenarios, -the method of 
this invention develops candidate strategies. To define 
a strategy, DMs select, out of large number n of 

35 variables x[ql [H] of the SSLP model, a small number K, K 
being substantially less than n, of key (or strategic) 
decision variables, say, first K variables xlq] [H] (q « 
1,..., K) . These K strategic variables are the same for 
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all scenarios and for all SSLP scenario models. They 
correspond to the most important decisions or to 
important immediate and irrevocable decisions. A 
candidate strategy I (I « 1,..., NST) is then defined as 
5 a vector of length K of ntmeri^al values for K strategic 
variables, with NST being the e^^ected nxamber of 
candidate strategies to be obtained from classification 
of the values of K strategic variables into segments* 

The centerpiece of the structure of this invention 
10 is the "outcome matrix" (OM) , which is a three- 
dimensional array of outcomes of risk- related activities, 
systematized in the "strategy versus scenario versus 
risk* related activity" format. The dimensions on the 
three corresponding axes are HST x S x where NST is 
15 the number of strategies, S is the number of scenarios, 
and R is the number of risk-related activities, so we 
must reserve in comiputer memory the necessary space for 
such an array, as well as for the regret matrix derived 
from the outcome matrix. (The number of dimensions may 
20 be greater than three.) 

The conpressed, highly aggregated results 
accumulated in the outcome and regret matrices serve two 
main purposes. First, they allow easy generation of non- 
predetermined, discretionary risk-limiting constraints 
25 and therefore help to form and modify candidate 

strategies, second, they allow the strategies to be 
analyzed and compared. 

If SSIjPb are used, this method has to generate 
strategies from S solutions of SSLP models. In 
30 coii5»rison with using an MSLP, this may be a drawback, 

because conc>romiee strategies resulting from MSLP tend to 
fit diverse conditions better. Also, the values of 
strategic variables in multiple solutions of an SSLP tend 
to differ somewhat, even if these solutions are close, so 
35 strategies have to be defined not by exsKrt values of 

variables, but by value segments where they belong. In 
using the method vrith an SSLP, DMs expect that S, the 
nuniber of scenarios, is sufficiently large, so at least 
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one of S BOlutlone will generate a good candidate 
strategy. 

The method of constructing candidate strategics and 
filling in the outcome matrix via SSLP models comprises 
5 three parts. First, it finds both eeasible and 

satisfactory solutions of models (3) . (The difference 
between ""feasible" and "satisfactory" is defined here as 
follows. A "feasible" solution meets all predetermined 
boundary limits on risk -related activities, which are 

10 identical for all scenarios. A "satisfactory" solution 
meets the predetermined limit, but also meets additional, 
"discretionary" limits on individual outcomes of specific 
risk- related activities under some specific scenarios, if 
DMs find these outcomes excessive or wish to find and 

15 explore a tradeoff.) Second, it splits the total value 
range for each of K strategic variables into one or more 
segments, assigns the values of these variables in each 
solution of an SSLP scenario model to the -relevant 
segments, and then forms multiple candidate strategies 

20 based on different combinations of these segments that 
have occurred in S solutions. Finally i the method 
adjusts the candidate strategies <o diverse conditions of 
all S scenarios and records the outcomes in the format of 
an outcome matrdLx. 

2S The first part of the method of using SSLP in 

forming the strategies for a physical system under 
uncertainty can easily be followed on the flow chart of 
FIG. 2. After generating the initial symbolic model (3) 
and storing it in the computer memory together with the 

30 values of the relevant parameters, predetermined boundary 
limits, aind value segments for strategic variables, Mis: 
form s scenarios - step 49; 

in turn, formulate the computational SSLP model (3) 
for each H (H - l,..., S) , which includes R risk-limiting 
35 constraints for all risk- related activities • st:ep 53; 

solve the model by any appropriate con^uter LP 
system - step 55; ' 
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record in the coniputer memory the valujss o£ outcomes 
obtained in the optimal solutions of model (3) £or each 
scenario H and £or each of R risk- related activities - 
step 57; 

5 if a feasible solution of (3) is not obtainable for 

at least one it means that some predetermined boundary 
limits are too tight, so modify the boundary limits and 
the corresponding rislc- limiting constraints and start 
solving again the resulting models (3) - steps ^7 and «9; 

10 continue the procedure of modification of boundary 

limits until feasible solutions of models (3) are found 
for all H - steps 51 to 59, 67 and ^9; 

decide whether the obtained set of out cranes is not 
yet fully satisfactory for all such activities under all 

15 conditions, that is, whether these outccoies, although all 
within the boundary limits, should not be further 
"discretionally" limited for some specific scenarios and 
some specific activities - steps 61, 63, and €5; 

if so, add the desirable "discretionary* lind.<8 and 

20 the corresponding risk- limiting constraints and solve 

again those of the resulting models (3} for the relevant 
scenarios and continue the procedure of adding or 
modifying the discretionary limits until satisfactory 
solutions of models (3) are found for all relevant 

25 scenarios - steps 53 to 61, 63 and ^5; and 

obtain solutions of models (3) with the final set of 
discretionary constraints for all S scenarios - steps 53 
to 61, 71 and 73. 

The second and third parts of the method of using 

30 SSLP in forming the strategies for a physical system 
under uncertainty can easily be followed on the flow 
chart of FIG. 3. After finding the values of K strategic 
variables in all S solutions of scenario models (3) and 
classifying these values in each solution by placing them 

35 in the value segments they belong to, DMs: 

define each set of segments for R variables 
encountered in at least one solution as a candidate 
strategy, forming NST (NST & l> candidate strategies. 
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where NST is already not an expected but Icnown number of 
strategies - steps 75, 77, and 79; 

add "strategy definition" constraints (3.€) to model 
(3), in turn forming models (3A) for each of the 
5 candidate strategies 1 (I « 1, ... , NST) and for each 
scenario with each of the models (3A) consisting of 
expressions (3.0) through (3.6): 

LBtll s XK[H] <UB[I], (3.6) 
where XK[HJ is a vector of strategic variables (where 
10 XK[H] is a subset of the set X[H]), while LB [I] and UB[I] 
are lower bounds and upper bounds for the relevant value 
segments, respectively, for strategy I (1=1,... , NST) 
* steps 83, 85, 87, 89, and 91; 

for each candidate strategy I, solve S scenario 
15 models (3A) - step 93; and 

if all S solutions are feasible, consider strategy I 
a viable candidate J and fill its out*cames for R rislc- 
related activities and for S scenarios into the 
corresponding columns of the ou-tccxne matrix - steps 95, 
20 105, 107, 109, and 111; 

otherwise, switch to analyzing the next candidate 
strategy (I + 1) - step 97. 

After con^leting the procedures displayed in FIG. 2 
and FIG. 3, the method with SSLP proceeds to the strategy 
25 selection procedure to be described later in connection 
with PIG. 12. 

To construct candidate strategies and to fill in the 
outcome matrix via an MSLP model, this invention uses a 
different method. The MSLP model generates not a 
30 multitude of solutions, but a single solution that is a 
conqpromise for diverse scenarios. However, similar to 
SSLP (3), from the very beginning, MSLP (4) also includes 
risk* limiting constraints for the predetermined boundary 
limits. 

35 The MSLP model (4) contains expressions (4.0) 

through (4.8) , %diere SN » S x n, SR « S x R, and 
SK » S X K: 
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10 



IS 



20 



25 



30 



35 



Find (SN X 1) vector X, X 
and (SR x 1) vectors 

XR, XR - (XR^, . 

W# W " * • • • 

Y, Y o lY^, ... 
and (SK x 1} vectors 

XK, XK c (XR]^. . 

XKP, XKP - (XKP]^ 



(X 



1' 



Ye), 



• • 9 

and XKM, XKM » (XXmj^* . 
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to A[H]X}| c. B[H), 
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1, -. 
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(4.2) 




Xj, * 0. 


<H - 
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. S) 


<4.3) 






(H - 


1, . . 


. S) 


(4.4) 




Yh . 0. 


(H - 






(4.5) 


XKh - ^ 


lOCj^ + XKPjj - XKMjj - 0, 


(H - 


2, . 


, . . S) 


(4.6) 


0 « XKpj] s UB, 


(H = 


2, . 


... S) 


(4.7) 




0 * XKotji s UB. 


(H - 






(4.8) 



In effect, the MSLP model (4) consists of S 



siibmodels (4.1) 
" Intersubmodel " 



- (4.5) for S scenarios H and of K(S - i) 
parts (4.6) - (4.8). The submodels are 
paired with the corresponding summands 28 (Hi of the 
objective function ZM described by (4.0). The submodels 
H correspond to SSLP (3) with the values of coii?>onents of 
their vectors C[H1 and BlH) and matrices A(H1 derived in 
the computer from the values of physical system 
parameters taken from the corresponding scenario H. 
Sxammands ZSlH] are weighted by coefficients (based on 
scenario probabilities) utHj of the corresponding 
scenarios H. 

In model (4) : X is an (SM x 1) vector of activity 
levels for n activities under S scenarios; XR is an (MR x 
1) vector of outcomes (where XR is a subset of the set X) 
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for risk- related activities, which are the last R -of M 
activities, xinder S scenarios; BL is an {NR x 1) vector 
of upper boundary limits for the same set of activities, 
invariant of scenarios; W and Y are, respectively, (NR x 
5 1) vectors of "slack variables* and "artificial 

variables" that transform risk^limiting inequalities into 
equations under S scenarios; and F is an (1 x NR) vector 
of fines f [ql (q - n - R ♦ 1, . . . , n) that are imposed 
on the excessive portions ytql tHl of out<:omes that exceed 

10 the boundary limits bl [q] (q » n - R 1, ... • n) • 
Equations (4.2) are NR risk- limiting inequalities 
transformed into equalities by adding slack variables w 
and subtracting artificial variables Y, while (4.3), 
(4.4) and (4.5) are the non - negat ivity condi t ions for 

15 variables XR, W and Y. The NK- dimensional vector XK 

(where XK is a subset of the set X) represents the values 
of K strategic variaibles under S scenarios, and equations 
(4.6) provide for the values of the K key variables being 
within close proximity of each other under all S 

20 scenarios H of the MSLP model (4) . Variables XKP and XKM 
characterize the "intersubmodel" differences between the 
values of strategic variables. Inequalities (4.7) and 
(4.8) provide for, on the one hand, the non-negativity of 
these variables and, on the other hand, the allowed upper 

25 bounds (UB) on their values. These bounds may differ 

between strategic variables, but are assumed here to be 
invariant of both of these variables and scenarios. 

FIG. 4 is a flow chart of a method, in accordance 
with the invention, for using MSLP to form strategies . 

30 The method begins with step 115, which stores an init:ial 
symbolic MSLP model (4) in computer memory, it continues 
with step 117, which enters parameter values for S 
scenarios. Next comes step 119, in which predetermined 
boundary limits are entered into memory, and step 121, in 

35 which scenarios are formed. Following that comes step 
123, which stores in memory a computational MSLP model 
(4) , smd step 125, which sets J equal to zero. In step 
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127, the MSliP model is solved, after which J is 
incremented by l in step 129 . 

After step 129 comes step 131, which defines 
candidate strategy J, and step 133, which records 
5 outcomes in computer memory. Prom there, the method 
moves to step 135, lAich compares outcomes with any 
predetermined and discretionary boundary limits, and 
decision point 137, which determines whether or not the 
strategy is completely acceptable. If the answer is yes, 
10 the method moves directly to final step 139, in which the 
selected strategy is in^lemented in the physical system 
or systems. If the answer is no, the method moves to 
decision point 141, which determines whether or not a 
sufficient number of new candidate strategies has been 
15 selected. If the answer is yes, the method moves to step 
151, which is the strategy selection process which will 
be described later and illustrated in FIG. 12. If the 
answer is no, the method moves to stap 143, which adds or 
changes discretionary constraints on outcomes, and from 
20 there to step 145, which solves the expanded model (4A) 
with discretionary constraints on outcomes. 

Following step 145, the method moves to decision 
point 147, which tests the feasibility of the model 
solution. If the answer is no, the method returns to 
25 step 143. If the answer is yes, the method moves to step 
149, which fills "strategy J" columns of an outcome 
matrix, and then back to step 129. 

The method of using MSLP for formation of strategies 
can be easily followed on the flow chart of PIG- 4. 
30 After forming initial syiribolic model (4), entering the 

values of parameters, and predetermined boundary limits, 
forming scenarios, and -forming and storing the 
computational MSLP model (4) in steps 115, 117, 119, 121, 
and 123, DMs: 

35 solve the MSLP model (4) by any appropriate computer 

LP system - step 127; 

define the vector of values of key variables in the 
optimal solution of (4) as Strategy J - step 131; 
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record in the computer memory the values o£ 
outcomes/ defined as OtH] tJl Ir] (H « 1, . ♦ . , S; r - i, . . • , 
R) , obtained in the optimal solution of Strategy J in 
submodels H of (4) for each of R risk- related activities 
5 r - step 133; 

conpare, for each (H, r) pair (H « S; r - 

1. . . . . R) , OtH] IJl trl with ref«rence8< targets and 
predetermined or discretionary limits, and deciding 
whether the obtained set of 0[H1 [J) [r] is fully 
10 acceptable for all H and r - step 135; 

if this set of 0[H] [Jl [rl , under all H and r {H - 
1, . . . , S; r = 1, . . . . R) , is both fully acceptable and 
superior to any other known strategy, the method stops 
the further developooent of candidate strategies and goes 
15 to implementation of strategy J - step 139; 

otherwise, DMs consider whether they have enough 
candidate strategies for their selection of a strategy - 
decision point 141; 

if they have, go to the "select strategy" process 
20 (see PIG. 12, which will be described later) and to 
strategy implementation - step 151; 

otherwise, transform in the computer . memory the MSLP 
model (4) into model (4A) , appending 'to it discretionary 
risk- limiting constraints and other modifications that 
25 restrict all values of 0(H1 [Jl trl to acceptable 

discretionary limits in5)OBed for individual scenarios and 
activities, or modifying such already existing 
constraints, to obtain a new strategy (J + i) with 
different outcomes 0[H] U+l] [rl • step 143; 
30 solve the MSLP model (4A) with discretionary risk- 

limiting constraints on the outcomes of risk- related 
activities by any appropriate ccanputer LP system - step 
145; 

check whether the resulting solution of MSLP (4A) 
35 with discretionary constraints generates a feasible 
candidate strategy • decision point 147; and 

if it does, fill the "Strategy J" columns of all 
two-dimensional •scenario vs. strategy" S x NST 
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subBiatrlces r (r « . . . , R) of OM wit:h 0(H] [J] [r] and 
proceed to constructing strategy (J ♦ 1) - step 149; 

otherwise « return to adding or changing 
discretionary constraints on outcomes - step 143. 
5 In spite of new^ efficient approaches to dealing 

with con^lexlty of decision-making that are incorporated 
in this invention, the nxiniber of scenarios or strategies 
still may turn out to be too large for solving the 
problem with available computer resources in ^he required 
10 time. In that case, each of these numbers can be reduced 
by grouping. There are two major forms of such grouping. 
One is "clustering," a well-Jcnown technique for grouping 
similar objects. The second is using the existing and 
known functional or statistical relationships between 

15 objects. 

Every method of clustering has its drawbacks. Using 
any specific method may lead to allocation of objects to 
the wrong clusters. In usual applications of clustering, 
such errors are difficult to find and eliminate. This 

20 invention is different: although it may perfectly well 
deal with object groups, it is especially interested in 
the "worst" scenarios and strategies. Therefore both 
objects with ambiguous cluster allocation and the worst 
objects of clusters can be separated and treated as 

25 individual entities - either from the very beginning or 
upon obtaining computational results. At some stages of 
the method, clusters or other groups that deserve special 
attention CcUi be corqpletely "unclustered, * and later 
clustered again, in old or new combinations . 

30 In PIG. 5, the method begins with step 153, in which 

an initial symbolic SSLP model (3) is stored in computer 
memory. Next comes step 155, in which parameter values 
for S scenarios are entered, and then step 157, in which 
predetermined boundary limits are entered. Then, step 

35 159 enters value segments for K strategic variables, and 
step 161 groups or regroups S individual scenarios into T 
typical scenarios. Prom there, step i€3 sets H equal to 
1. 
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Next in FIG. 5 comes step 165, which generates in 
computer memory a computational model for scenario H. 
Prom there, step 167 solves the ccsnputatibnal model for 
scenario H and step 169 obtains and records outcome 
5 levels in memory* Following step 169 comes decision 

point 171, irtiich determines whether or not the recorded 
solution is feasible. If the answer is yes, decision 
point 173 determines whether all outcomes are 
satisfactory. If the answer is no, the method moves to 

10 step 175, which adds or modifies discretionary 

constraints limiting imsatisfactory outcomes, to step 
177, which sets H equal to 1, and back to step 165 for 
further computational model generation. 

If the answer at decision point 173 in FIG. 5 is 

15 yes, the method moves to decision point 179, which 

determines whether or not H is equal to T. If the answer 
is no, the method moves to step 181, which increments H 
by 1 and returns to computational model generation step 
165. If the answer at decision point 171 is no, the 

20 method moves to step 183, which modifies at least one of 
the predetermined boundary limits, to step 185, which 
modifies the risk-limiting constraints, and back to step 
163, which sets H to 1. Finally, if the ansifer at 
decision point 179 is yes, the method moves to decision 

25 point 187, which determines whether regrouping is 

desired. If the answer is yes, the method returns to 
step 161. If the answer is no, the method goes to step 
189, which returns to step 75 in 3. 

The method for finding the solution for a physical 

30 system under uncertainty via the use of SSLP models with 
clustering can easily be followed on the flow chart of 
FIG. 5. It can be seen that the flow chart in PIG. 5 is 
almost identical to the flow chart in FIG. 2. The 
exception is that the "FORM S SCENARIOS" stage is 

35 replaced by the "GROUP OR REGROUP * INDIVIDUAL* INTO 

'TYPICAL* SCENARIOS" (Step 161), SO that one method loop 
returns to that operation. 
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Another approach to solving decleion-malcing problems 
where the nuinber of scenarios is too large to -cover with 
a single MSLP model (4) is by solving many MSLP models 
(4> of smaller size. Let the total number of scenarios 
5 equal S and let the largest MSLP model (4) that can be 
solved in acceptable time have WS scenarios. For 
slnplicity^ let also NDMS - S/WS be an integer. The 
method of solving the problem by this technique can be 
followed easily on the flow chart of FIG. 6. 

10 In FIG. 6, the method begins with step 191, in whi<:h 

an initial synibolic MSLP model (4) is stored in computer 
memory. Next, the method moves through step 193; in 
which parameter values for S scenarios are -enter^ed into 
memory, step 195, in which predetermined boundary limits 

15 are entered into memory, and step 197, in whi<:h S 
scenarios are formed. Next, step 199 clusters S 
scenarios into T groups of WS scenarios. From there, 
step 201 sets J equal to zero, where J is as defined 
heretofore. 

20 Prom step 201 in FIG- 6, the method moves to step 

203, setting H equal to 1, where H is an index of ^he 
scenario group being tested, and to step 205, which 
generates a computational MSLP model (4) or (4A) <£or 
group H. After step 205 comes step 207, which solves the 

25 model (4) or (4A) for group H, and step 209, which 

obtains and records outcome levels. After step 209 -^omes 
decision point 211, which determines whether or not the 
solution is feasible. If the answer is yes, the method 
moves to decision point 213, which determines whether all 

30 outcomes are satisfactory. If the ainswer there is yes, 
the method moves to decision point 215, whi^h determines 
whether or not H is equal to T. If the answer at 
decision point 215 is no, the method moves to step 217, 
which increments H by 1, and proceeds back to model 

35 generation step 205. If the answer at decision point 211 
is no, the method moves to step 221 to modify the 
predetermined boundary limits and to step 222 to modify 
the risk- limiting constraints before returning -to step 
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203, setting H equal to 1. If the answer at decision 
point 213 is no, the method goes to step 219, to add or 
modify "discretionary" constraints limiting 
unsatisfactory outcomes, to step 220, setting H equal to 
5 1, and back to step 205. 

In FIG. 6, if the answer at decision point: 215 is 
yes, the method moves through step 216 to FIG. 6 A, which 
begins with step 218 and moves to step 223, choosing the 
WS worst scenarios. Also necessary for conpletion of the 

10 method illustrated in FIG. 6 are point A, which comes 
from step 227 in PIG. 6A and point B, which comes from 
step 229 in FIG. 6A. 

From FIG. 6, a yes answer at decision point 215 
connects to step 223 in FIG. SA, in which the WS worst 

15 scenarios are selected. From there in FIG. 6A^ -the 

method moves to step 224, forming an fGLP model (4) for 
WS scenarios. Next, step 225 solves the model and moves 
to decision point 227 to determine whether or not the 
solution is feasible. If the answer is no, the method 

20 moves to point A and back to step 221 in FIG. 6. If the 
answer is yes, the method moves -to decision point 229 to 
determine whether or not all outcomes are satisfactory. 
If the answer is no, the method moves to point B and back 
to step 219 in FIG. 6. If the answer is yes, the method 

25 moves to step 231 to obtain the values of X strategic 
variables. 

Prom step 231, the method in FXG. 6A moves to step 
233, where J is incremented by one. From there, it moves 
to step 239 to define strategy J and formulate "strategy 

30 definition" constraints for strategy J. Prom step 239, 
the method moves to step 241, where H is incremented by 
1, to step 243, %irhere a model (4) for group H is added to 
the strategy definition constraints for strategy J, and 
to step 245, where the model (4B) for group H is solved. 

35 From step 245, the method moves to decision point 247 to 
test whether or not the solution is feasible. Xt the 
answer is yes, the method moves to decision point 249 to 
determine whether or not H is equal to T. If the answer 
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at decision point 249 is yes* the method moves 4:o point C 
and thence to PIG. 6B in a manner which will be described 
later. If the answer at decision point 249 is no, the 
method moves back to step 241 in PIG. 6A. 
5 If the answer at decision point 247 is no, the 

method in FIG. 6A moves to step 251, where J is 
decremented by 1. Prom step 251, the method moves to 
step 253, which modifies predetermined or •discretionary" 
restraints on outcomes, and to step 25S, which sets H 

10 equal to 1. From step 255, the method moves back to step 
224, to form an expanded MSLP model for WS scenarios. 

In FIG. 6B, step 257 is driven from point C, which 
comes from a yes answer at decision point 249 in FIG* 6A. 
After step 257, step 259 fills the "Strategy J" columns 

15 of an outcome matrix. From there, the method moves to 

decision point 261 to determine whether or not to ^all a 
halt to the generation of new candidate strategies. A 
yes answer leads to step 263, which is ^he "select 
strategy" process yet to be described in connection with 

20 FIG. 12. Following step 263, the method comes to an «nd 
with step 265/ which implements the selected strategy in 
the physical system or systems. 

A no answer at decision point 261 in FIG. 6B causes 
the method to move to step 267, which changes the 

25 con^osition of the group of WS scenarios. From there, 
the method moves to point D, which goes to step 255 in 
PIG. 6A. 

Similarly to the method above that uses an SSliP 
model (3>, the method currently being described consists 

30 of three parts. First, it clusters S scenarios into NUMS 
groups of WS scenarios each, forms model (4) for -each 
group, finds a solution of that model (4) that is both 
feasible and acceptable, and obtains axui records the 
values of outcomes of all risk- related activities. 

35 Second, the method selects a group of WS scenarios that 
are the worst in some outcomes, and repeats the above 
procedure for that "worst group" - that is, forms model 
(4) for the group, finds a solution of that model (4) 
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that is both feasible and acceptable « and obtains and 
records the values of outcomes of all risk- related 
activities. The method also obtains the values of 
strategic variables in the solution of model (4> and 
5 defines the corresponding candidate strategy. Finally, 
the method adjusts the candidate strategy to diverse 
conditions of HUMS scenario groups (t:hat is, of all S 
scenarios) and records the outcomes in the format of an 
outcome matrix. This process can be repeated t.o form as 
10 many candidate strategies as desired, starting from 

different scenario groups and different discretionary 
cozwtraints. 

In the first part, as illustrated in PIG. €, at<er 
forming the initial symbolic model (4), entering the 
15 values of parameters, predetezmined boundary limits, and 
value segments for strategic variables , DMs : 

form S scenarios - step 197; 

split S scenarios into T groups, ^each conprised of 
WS scenarios - step 199; 
20 for each of these groups, in turn, generate and 

store computational NSLP model (4) and solve model (4) by 
any appropriate computer LP system - steps 205 and 207; 

if any of T solutions of confutation^.! MSLP models 
(4) is infeasible, modify the values of predetermined 
25 boundary limits and the corresponding risk* limiting 

constraints, and solve all T models (4) again - steps 
211, 221, 222, 205, and 207; 

obtain and record the values of outcomes - step 209; 

and 

30 if, for some scenarios and some activities, ^he 

outcomes are not considered satisfactory, isqpose or 
modify (for the scenario groups and activities involved) 
additional "discretionary" risk- limiting <:on8t^aints and 
solve computational models (4A) again for the relevant 

35 scenario groups - steps 219, 220, 205, and 207. 

In the second part, as illustrated in PK3S. ^ and 
6A, when all outcomes for all scenario groups are 
satisfactory, DMs: 
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form a new scenario group tram W3 scenarios where 
the ou teenies are the worst, form an MSLP model (4) €or 
that group and solve the model by any apprc^riate 
coniputer LP system 1 - steps 223, 224, and 225; 
5 again, <:heck whether the solution is feasible and. 

if not, modify the predetermined limits - steps 227, 221, 
and 222; 

if desirable, modify the discretionary constraints 
and obtain both feasible and satisfactory solution - 
10 steps 229, 219, 205, and 207; 

obtain and record the values of outcomes, including 
the values of strategic variables • steps 209 and 231; 
and 

formulate "strategy definition" constraints - step 

15 239. 

in this part of the method, screening out of bad 
solutions is accomplished by a combination of three types 
of constraints: constraints based upon predetermined 
boundary limits, discretionary constaraints , and "strategy 

20 definition" constraints that limit the values of 
strategic variables. 

in the third part, model (4) f or jeach of T groups of 
WS scenarios is, in turn, added to "strategy definition" 
constraints to form model (4B) for that group. 

25 Conputational models (4B) are solved in turn for all T 
groups. The values of outcomes in risk- related 
activities for all S scenarios are obtained from these 
solutions and fill the coltimns for the current strategy 
in the outcome matrix. The method starts a new candidate 

30 strategy by changing either the discretionary constraints 
or the "worst scenario" group. 

As indicated above, this method is a complex part of 
the invention, but it is also the best for dealing with 
very large and complicated problems. Even if the size of 

35 a problem is so large that no more than two scenarios can 
be put into an MSLP. its capacity to find compromise 
solutions, along with its ability 4:o form many 
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combinations of scenarios « permits a sharp inprovement in 
the quality of strategies. 

This concludes the section of the "Detailed 
Description" about aspects of the invention related to 
5 f onnulating and solving LP models in constructing 

candidate strategies. The next part refers to methods of 
selecting a strategy to be implemented in the physi<3l 
system. 

Most state-of-the-art methods of selecting a 

10 strategy on the basis of a payoff (outcome) matrix are 
described in the previously quoted "Barnes and Solutions 
Chapter 13. This invention offers three new methods 
described below. State-of-the-art coiq3uter methods refer 
to two-dimensional payoff matrices that contain "strategy 

15 versus scenario" outcomes. Two-dimensional matrices are 
sufficient if there is only one risk-r-elated activity, 
such as the financial results of a project. 

FIG. 7 demonstrates such a two-dimensional payoff 
matrix for a four- strategy, six-scenario problem. The 

20 single risk- related activity is here the cost of a supply 
chain production and distribution plan, which is 
developed by using MSLP models (4) and <4A) in a real 
world problem. The matrix entries are the costs of the 
strategy under different "strategy versus scenario" 

25 conibinations (in dollara) . 

In addition to the true payoff matrix shoim in the 
first six lines and four columns of cose numbers, PiG. 7 
contains four aggregate cost results (BEST -COST, WORST 
CX)ST, EXPECTED COST, and EXPECTED IR COST) for each 

30 strategy., and one aggregate (BEST SCENARIO <X>ST) for .^ach 
scenario. IR stands here for "insufficient reason"; this 
term is explained below. W,^ECTED COST is computed by 
Eq. (5) and EXPECTED IR COST by Eq. (6) : " 
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S 

EXPECTED COST (J) " T u[H]0[H](J]; (5) 

H»l 

S 

EXPECTED IR COST (J) = (E 0W1[J]) / S; (6) 

where 0[HI [J] la the cost payof£ Cor strategy J under 
scenario utHl is the probability of scenario H, and S 
is the number of scenarios. 

In payoff -based criteria, four strategy cost results 
are used for comparing the strategies » screening out the 
worst strategies, and, finally, selecting the best one. 
The best scenario results are used for computing the 
regret matrix. 

For amy "strategy versus scenario" combination, 
regret is computed by Bq. (7) : 

RBGRBT[H11J1 - BESTIHJ - 01HJ[J1, (7) 
where: BEST(H} is the BEST SCEMfkiao COST for scenario H 
and O IHl ( Jl is the payoff for strategy J under scenario 

For payoffs of FIG. 7, the regret matrix is shown in 
FIG. 8. In each scenario row, at least one entry (for 
the strategy that is the best under that scenario) equals 
zero. For any other strategy, REGRET is the amount of 
opportunity lost because the strategy, if selected, turns 
out not to be the best under scenario H. Again, in 
addition to the true regret matrix in the first six lines 
and four columns of regret numbers, FI<3. 8 <:ontain8 four 
aggregate regret results for each strategy. EXPECTED 
REGRET is computed by Bq. (8) and EXPECTED IR REGRET by 
Eq. (9) : 

S 

EXPECTED REGRET (J) = S u [H| REGRET [HH J] ; (8) 

H»l 
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EXPECTED IR REGRET (J) » (E REGRET (HI [J] ) / S (9) 

5 where REGRET IH] [J] is regret for strategy J under 

scenario H, u(H] is the probability o£ scenario H, and S 
is the number of scenarios. 

Three new regret -based methods ^o be introduced in 
this invention are: 
XO (a) OPR - The Optimism- Pessimism Index (Regret) 

Criterion; 

(b) PIR - The Partial Ignorance (Regret) Criterion; 

(c) MIRR * The Modified Insufficient Reason <Regret) 
Criterion. 

15 The previously known payoff -based counterparts for 

t:hese methods are: 

For OPR - OP, or The Optimism- Pessimism Index 

Criterion; 

For PIR - PI, or The Partial Ignorance Criterion; 
20 For MIRR - MIR, or The Modified Insufficient Reason 

Criterion* 

For all three methods and their known counterparts, 
the general methodology has many common points. It is as 
follows. 

25 If decision-makers (DMs) are very pessimistic, «iey 

can centre the strategies just on the basis of ^helr 
worst results, which are WORST <X}ST for the payoff 
methods and WORST REGRET for the regret methods. This 
corresponds to previously kno%m "minimax payoff" and 

30 "minimax regret" criteria, r^espectivaly. In other words, 
DMs choose the strategy that has the best WORST and thus 
try to protect themselves against disastrous outcomes. 

If DMs are very optimistic, they can conrpare the 
strategies just on the basis of their best results, which 

35 are BEST COST for the payoff methods and BEST REGRET for 
the regret methods. This corresponds t.o previously known 
"maximsuc payoff" or "maximax regret" criteria. 
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respectively. In other words, DMs choose ^he strategy 
that has the best BEST. 

Decision Science offers the previously loioim 
synthetic OP (•optimism-pessimism index") criterion to 
5 retain the flexibility of DMs* and their ability ^o 

change the selection criteria at will, depending on their 
judgment of the conditions of a specif i<: problem. The OP 
criterion is applicable to payoffs only. The 
corresponding new method of this invention, OPR« is 
10 applicable to regrets only. 

Both criteria consider as unreliable any information 
on scenario probabilities, and therefore compare 
strategies J only on the basis of BEST (OOST or REGRET, 
respectively) , WORST (COST or REGRET, respectively) , and 
15 an "optimism-pessimism index" ALPHA. ALPHA incorporates 
both objective Icnowledge, including the existing 
forecasts 4 and subjective judgment of DMs. 

The first new method involves the OPR 
criterion, which is: 
20 Select J that minimizes 

ZOPRCJ, ALPHA) « (1 - ALPHA) ♦ BEST REGRET (J) ♦ 
ALPHA ♦ WORST REGRET (J) , 
0 s ALPHA s 1 (J»l«..., NST), (10) 
Where the value of ALPHA is specified by the user and HST 
25 is the nuskber of strategies. ZOPR (ALPHA) is here the 
expected regret for strategy J if probability ALPHA ia 
assigned to the scenario that is the worst for J, 
probability 

(1 - ALPHA) to the best scenario, and zero probabilities 
30 to all other scenarios. The previously known OP 

criterion is similar, with ZOPR replaced by ZOP and 

REGRET replaced by COST. 

The extreme value of ALPHA = 1 corresponds to the 

pessimistic outloo)c of minimax regret, while ALPHA « 0 
35 corresponds to the optimistic outloo)c of maximax. Since 

an arbitrary choice of a specific value of ALPHA does not 

inspire confidence, a set of different values of ALPHA 

can be used. The values of ZOPR (J. ALPHA) under varying 
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values of ALPHA/ as they are described by (10), can be 
represented by a straight -line segment on '^he (0,1) 
interval; see FIG. 9. 

FIG. 9 is a decision graph illustrating the 
5 optimism-pessimism index (regret) for strategies B« C, 
and D using the OPR criterion. The ordinate is regret, 
measured in dollars, and the abscissa is ALPHA, measured 
from zero to unity. Strategy A is omitted, since it is 
obviously much worse. In PIG. 9, the aOPR line for 
10 strategy B is: 

ZOPR(B, ALPHA) « (1 • ALPHA) * 0 4- ALPHA * 10,008. 
The ZOPR line for strategy <: is: 

ZOPR(C, ALPHA) - (1 - ALPHA) * 4,9S0 + ALPHA * 6,541. 
These two lines intersect at the value of ALPHA 

15 where the ZOPRs for both strategies are jequal, or at 
ALPHA « 0.5681. This means that if DMs consider the 
probability of the best sceziarios (which, for both 
strategy B and strategy happen to be scenarios 1 and 
5) as less than 0.5881, they should select strategy B 

20 because it has, up to that value of ALPHA, smaller values 
of ZOPR. On the other hand, if they are pessimistic and 
consider ALPHA closer to i, they should select strategy 
C. 

FIG. 9 presents the "OPR efficient frontier'' for 
25 this example. The efficient frontier is the set of best 
strategies on the (0, 1] range. It comprises two bold 
straight-line segments: from (0, 0) to (0.5881, 5,886) 
and from (0.5881, 5,886) to (1.0, 6,541). Bach strategy 
on the efficient frontier is characterized by a 
30 Bubinterval of the [0, 1] range where this strategy is 

the optimal choice. The efficient frontier incorporates 
valuable information about the relative merits and 
demerits of any strategy: 

• The composition of the subset of strategies that 
35 form the efficient frontier. 

• The width of the subinterval supporting each 
optimal strategy. 
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• The order o£ optimal strategies from the 
optimistic end of the [0, 1] interval to the pessimistic 
end. 

• The difference between the regrets for the best 
5 strategy and other strategies ^ which shows the possible 

impairment of results in choosing a non- optimal strategy. 

The OPR efficient frontier consists of two 
strategies, B and C. The "optimistic part" of the 10, 1] 
interval, where strategy B is optimal, has a width of 
10 0.5881, while strategy C is optimal on a 0.4119 

subinterval. At ALPHA - 0.5, the difference between the 
regrets of C and B equals 742, while at ALPHA » 0.7 it 
equals -942. 

The second new method involves the PlR criterion, 
15 very similar to the OPR criterion. The PIR criterion is, 
however, more cautious: while leaving intact the 
"pessimistic end" - that is, WORST REGRET (J) - of the 
straight* line segment that represents strategy J, it 
replaces the •optimistic end" by the EXPBCTB> REGRET [J] . 
20 The PIR criterion is: 

Select J that minimizes 
ZPIR(J, ALPHA) - (1 - ALPHA) * EXPECTED REGRET<J) ♦ 
ALPHA * WORST REGRET (J) , 
0 s ALPHA s 1 (J-1,..., NST), (11) 
25 where the value of ALPHA is again specified by the user. 

(The previously )cnown PI criterion is similar, with ZPIR 
replaced by ZPI and REGRET replaced by COST.) 

The extreme value of ALPHA •» l again corresponds to 
the pessimistic outlook of minimax regret, while ALPHA « 
30 0 corresponds to a more optimistic outlook of expected 
regret . In other %fords , DMs do not want to be -carried 
away by overly rosy expectations of potential BEST 
REGRET; they therefore lower their expectations and 
assume that the best they <:an hope for is EXPECTED 
35 REGRET. At the same time, the real regret can be much 
lower, do%m to WORST REGRET. The more reliable the 
forecast scenario probabilities, the closer DMs can move 
toward the slitple expected regret criterion by decreasing 
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the value of MiPHA. In contrast, the less DMs trust the 
information available, the more they should be inclined 
to be conservative. At ALPHA « l, the PIR method is 
again reduced to the pessimistic minimax regret 
5 criterion. 

The values of ZPIR(J« ALPHA) under varying values of 
ALPHA, as they are described by Bq.(ll), can be 
represented by a straight- line segment on the (0,1) 
interval; see FIG. 10, which is a graph of PIR efficient 
10 frontier. FX6. 10 is a decision graph illustrating the 

partial ignorance index (regret) for strategies B, C, and 
D using the PIR criterion. The ordinate is regret, 
measured in dollars, and the abscissa is ALPHA, measured 
from zero to unity. 
15 FIG. 10 also depicts only strategies B, C, and D. 

In FIG. 10, the ZPIR line for strategy B is: 

ZP1R(B, ALPHA) - (1 • ALPHA) * 2,742 
^ ALPHA * 10,008. 
The ZPIR line for strategy C is: 
20 ZPIR{C, ALPHA) - (1 - ALPHA) ♦ 5,754 + ALPHA ♦ €,541. 

These two lines intersect at the value of ALPHA » 
0.4649 where ZPIRs for both strategies are equal. If the 
probability of the best scenarios is less than 0.4649, 
strategy B should be selected, otherwise - strategy C. 
25 Since the PIR criterion is more realistic than OPR, the 
area of an optimistic choice of a strategy becomes more 
narrow. 

The third method involves the MIRR criterion used 
when the information on scenario pr<d>abilities is either 

30 absent or quite unreliable. In such cases, there is no 
sufficient reason to consider that the scenarios are not 
equally likely to occur, so the probability of jeach 
scenario is assumed to be 1/S, where S is . the number of 
scenarios. With BXPBCTRD IR RBGRBT computed by Bq. (9) , 

35 the MIRR criterion is quite similar to the PIR criterion. 
The MIRR criterion is: 

Select J that minimizes 
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ZIRR(J, ALPHA) - (1 - ALPHA) * BXPBCTED ZR KSG^STIJ) ^ 

ALPHA * WORST RBGRETtJ) » 
0 s ALPHA ^ 1 (iJ»l,..-, NST) . (12) 
The straight-line ZIRR segments based on Sq. (12) are 
5 8ho%m In FIG. 11 « which is a graph o£ MIRR efficient 
frontier. This time, the intersection of lines for 
strategy b and strategy C is at 0.3552. PIG. 11 is a 
decision graph illustrating modified insufficient reason 
(regret) for strategies C, and D using the MIRR 

10 <:riterion. The ordinate is regret, measured in dollars, 
and the abscissa is ALPHA, measured from zero to unity. 

The three new regret -based methods of this invention 
allow much more sopliisticated and sharp analysis in 
conparing and selecting strategies. This is demonstrated 

15 in the example depicted in FIGS. 7 through 11. The cost 
matrix in FIG. 7 has been generated in an atten^>t to 
reduce the rislclness of evolving candidate strategies - 
first strategy A, then B, then C. Strategy A is very 
risky. Compromise strategy B is generated by using model 

20 (4) and the corresponding multiscenario model. The next 
step is intended to reduce the high costs encountered 
under scenarios 2 and 4 - even at the pri^e of increasing 
costs under other scenarios. Imposing the needed 
discretionary constraints achieves that goal and obtains 

25 strategy C. Repeating this procedure with tighter 

discretionary constraints further reduces the costs under 
the two worst scenarios. It therefore seems that 
strategy D is less rislcy than strategy C, which in turn 
is less risky than strategy B. The cost results more or 

30 less confirm that impression: in moving from strategy C 
to strategy D, worst "costs decrease by 5,266, while the 
expected costs rise only 3,347 (see FIG. 7). 

Analysis of the regret results and of the regret - 
based efficient frontiers of FIGS. 9 through 11, however, 

35 brings one to a conpletely different conclusion. The 
regret outcomes for strategy D xinder scenarios 2 and 4 
indeed become the best. However, <:he price paid is too 
steep: regrets under the other four scenarios exceed 
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10,400 and under two of them -come to aroimd 12,000* 
Correspondingly, as can be seen, £or Instance, in FIG* 
10, strategy D is much worse than either strategy B or 
strategy C. Such incisive analysis is not possible when 
5 only previously known DS methods are used. The regret - 
based criteria are therefore invaluable in determining 
desirable limits on tightening discretionary constraints. 

To sinqplify and possibly improve the use of outcome 
and regret matrices in decision-making, the following 

10 difficulties (in addition to introducing new and iinproved 
DS criteria) have to be overcome: 

First, if, in addition to the outcomes of a first 
risk- related activity, such as cost, the problem has to 
deal with the outcomes of a second risk- related activity, 

15 such as the output of a dangerous by-product, a second 

two-dimensional outcome matrix results. The method thus 
generates a three-dimensional array of two matrices. The 
outcome matrix nay also have more than three dimensions - 
for instance, if the outcomes are classified by <the 

20 causes of uncertainty. 

These multidimensional matrices have ^o be 
transformed into two-dimensional ones where the state-of- 
the-art methods of Decision Science «an be aK>lied. 

Second, decisions are often made not by a single 

25 decision-maker, but by a group whose members may have 
different priorities and criteria for selecting a 
strategy. 

Third, some tables of the outcome matrix may be 
filled in not with numerical values, but rather wi^h 

30 qualitative characteristics such as "good" or 
•convenient." 

Fourth, there are several state-of-the-art and new 
criteria of Decision Science for selecting a strategy. 
They may give conflicting answers. 

35 This invention provides techniques for dealing with 

all these and similar difficulties. The key is in 
weighting the matrix entries. For instance, in 
converting a three-dimensional matrix into a two- 
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dimensional one, one risk-related activity, such as the 
financial results, is given priority as the xnain 
activity. The outcomea in all other risk*- related 
activities are converted into the additional outcomes of 
the main activity by multiplying them by conversion 
coefficients, or weights. The values of these 
coefficients are not likely to be known with accuracy. 
Consequently, several different values are used 
simultaneously for each conversion coefficient, so the 
number of scenarios considered in ^he matrix increases 
proportionately. 

The strategy selection method illustrated by the 
flow chart in PIG. 12 begins with step 2«9, where the 
purpose of the steps that follow is to process an outcome 
matrix and exclude strategies with unsatisfactory 
outcomes in any activity. From step 269, the method 
moves to decision point 271, where it is determined 
whether NUMDIM (the number of dimensions of the matrix) 
is greater than 2. If the answer is no, the method 
proceeds to decision point 273, where it is determined 
whether NUMGR (the number of members of the decision 
making group) is greater than 1. If the answer there is 
no, the method proceeds to step 275, where it is 
determined whether NUMQ (the number of types of 
qualitative characteristics that have to be transfocmed 
into quantitative values) is greater than zero. 

If the answer at decision point 271 in FIG. 12 is 
yes, the method moves to step 289 to specify weights WD 
and to step 291 to expand the matrix to matrix-D. From 
step 291, the method moves to decision point 273. If the 
answer at decision point 273 in FIG. 12 is yes, the 
method moves to step 285, where weights WG are specified, 
and to step 287^ where the matrix is expanded to matrix- 
G. 

If the answer at decision point 275 in FIG. 12 is 
no, the method moves to step 277 to specify weights WC 
and to step 279 to expand the matrix to matrix-C. If the 
answer at decision point 275 is yes, the method moves to 
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Step 293 to specify weights WQ and to step 295 to expand 
the matrix to inatrix*Q. 

Finally in PIG. 12, the method moves from step 279 
to step 281 and applies Decision Science criteria to 
5 exclude the worst strategies and to step 263 to use risk 
attitude and the subjective considerations of mfls, 
together with Decision Science criteria, to select the 
right strategy. 

The method of processing the outcome matrix or -the 
10 criteria for selecting a strategy by weighting can easily 
be followed on the flow chart of PIG. 12. After studying 
the payoff matrix and excluding those strategies that 
have unsatisfactory outcomes in any scenario or any 
activity, DMs: 

15 determine NUMDIM, the number of dimensions of the 

matrix; NUMSR, the number of members of the decision- 
making group; and NUMQ, the number of types of 
qualitative characteristics that have to be transformed 
into quantitative values - steps 271, 273, and 27S; 

20 if NUMDIM > 2, NUMSR > 1, and NUMQ > 0, develc^ a 

set of conversion weights for each of the respective 
transformations of the outcome matrix, as well as weights 
for different criteria for selection of a strategy - 
steps 285, 289, 293, and 277; 

25 use the weights to expand the outcome matrix as 

desired (MATRIX-D is the initial matrix after its 
transformation into a two-dimensional matrix; MA1^IX-<S is 
the matrix after group member differences are ^aken into 
account, MATRIX-Q takes into account qualitative result:s; 

30 and MATRIX- C smooths over the differences between 

different DS criteria) - steps 291, 287, 295, and 279; 
and 

apply the criteria of Decision Science, as well as 
considerations of risk attitude and sxabjective 
35 preferences of DMs, to select a strategy on the basis of 
both the original outcome matrix and its expanded version 
- steps 281 and 283. 
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It is to be understood that the enibodiinents of the 
invention which have been described are illustrative. 
KumerouB other arrangements and modifications may be 
readily devised by those skilled in the art without 
5 departing from the spirit and scc^e of the invention. 
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Claima: 

1. A conqputer method for mazxaging risk In multiple 
parameter physical Byetems performing interrelated 
activities, at least one of such activities being 

5 designated as a risk-related activity and having an 

outcome level which may fall outside of boundary limits « 
said method establishing a course of action (strategy) 
for said physical systems preventing any outcome levels 
for said risk- related activities from falling outside of 

10 said boundary limits and comprising the steps of : 

CD finding a set of satisfiable boundary 
limits (11-15) and developing (17-21, 25-29, 37) in 
conputer memory a multitude of feasible candidate 
strategies that satisfy said boundary limits; 

15 (2) finding the values of outccwies (23, 35) of 

said risk-related activities for each said feasible 
candidate strategy under relevant conditions by storing 
(17) and solving (19) a model of said physical systems in 
said computer memory; 

20 (3) identifying cazididate strategies from said 

outcomes of said risk- related activities and recording 
said outcomes in said con^uter memory; and 

(4) jointly applying in said computer memory 
multiple optimization criteria to said outcomes of said 

25 candidate strategies to aid decision-mEUcers in selecting 
an intplementable strategy. 

2. A computer method for managing risk in multiple 
parameter physical systems performing interrelated 
activities, at least one of such activities being 

30 designated as a risk- related activity and having an 

outcome level which may fall outside of boundary limits, 
said method establishing a course of action (strategy) 
for said physical systems preventing any outcome levels 
for said risk- related activities from falling outside of 

35 said boundary limits and conqorising the steps of: 

(1) finding a set of satisfiable boundary 
limits (11-15) and developing (17-21, 25-29, 37) in 
conputer memory a multitude of feasible cand id ate 
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Strategies that satisfy said boiindary limits by 
conqputational methods that use risk -limiting constraints 
(31-33, 39) ; 

(2) finding the values (23-35) of outcomes of 
5 said risk- related activities for each feasible <:andidate 

strategy under relevant conditions by storing (17) and 
solving (19) a model of said physical systems in said 
conqputer memory; 

(3) identifying said candidate strategies from 
10 said outcomes of said risk- related activities and 

recording said outcomes in said •computer memory; 

(4) jointly applying multiple optimization 
criteria in said coxnputer memory to outcomes of all 
candidate strategies to exclude candidate strategies 

15 which are inferior to others under said criteria; and 

(5) jointly applying multiple optimization 
criteria in said computer memory to outcomes of remaining 
candidate strategies to aid decision-makers in selecting 
an implementable strategy. 

20 3. A computer method for managing risk in multiple 

parameter physical systems performing interrelated 
activities, at least one of such activities being 
designated as a risk- related activity and having an 
outcome level which may fall outside of boundary limits, 

25 said method developing a course of action (strategy) for 
said physical systems preventing any outcome levels for 
said risk- related activities from falling outside of said 
boundary limits and conprising the steps of: 

(1) storing an initial model (11-17) of said 
30 systems in computer memory, said model including values 

for parameters of each of said physicail systems; 

(2) producing a solution (19) of said initial 
model in said computer memory; 

(3) deriving from said solution of said initial 
35 model a set of outcome levels (21-23) for said risk- 
related activities; 

(4) cosq[>aring said set of outcome levels (25- 
29) for said risk- related activities with said boxindary 
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limits and, if any of said boundary limits are violated, 
then deriving at least one risk* limiting derivative 
constraint from such conqparison; 

(5) expanding said initial model in said 
5 computer memory by adding said risk- limiting derivative 
constraints (31) to said initial model in said computer 
memory; and 

<6) producing a solution (33) of said expanded 
model with said added risk- limiting derivative 
10 constraints in said computer memory. 

4 . The computer risk management method of claim 3 
where said risk- limiting derivative ^oxistraints are in 
the form of mathematical equations and inequalities. 

5. The computer risk management method of claim 3 
15 where said risk- limiting derivative constraints are in 

the form of fines and penalties. 

6. A computer method for managing risk in multiple 
parameter physical systems performing interrelated 
activities, at least one of such activities being 

20 designated as a risk-related activity and having an 

outcome level which may fall outside of boundary limits, 
said method developing a course of action (strategy) for 
said physical systems preventing any of said outcome 
levels for said risk- related activities from falling 

25 outside of said boundary limits and cosnprising the steps 
of: 

(1) storing an initial model (11-17) of said 
systems in said computer memory, said model including 
both a set of decision variables and values for 

30 parameters of each of said physical systems, where 
strategic variables are a subset of said decision 
variables and where a strategy consists of a set of 
values specified for all said strategic variables; 

(2) producing a solution (19) of said initial 
35 model in said computer memory; 

(3) deriving from said solution of said initial 
model a set of outcome levels (21-23) for said risk- 
related activities; 
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(4) comparing said aet of outcome levels <25- 
29) for said risk-related activities with said boundary 
limits and, if the outcome level of at least one said 
risk* related activity falls outside said boundary limits 

5 for that activity, then deriving at least one risk- 
limiting derivative constraint from such comparison; 

(5) esqpanding said initial model in said 
computer memory by adding said risk- limiting derivative 
constraints (31) to said initial model in said computer 

10 memory; 

(6) producing a solution (33) of said expanded 
model with said added risk- limiting derivative 
constraints in said computer memory; and 

(7) identifying a strategy for said physical 
15 systems from the values of said strategic variables in 

said solution of said expanded model. 

7. The computer risk mauoagement method of claim 6 
in which said computational model of said physical 
systems is used to locate regions and groups of risky 

20 decisions and strategies euid to cut off said regions and 
groups by adding to said model derivative constraints. 

8« A computer method for managing risk under 
certainty in multiple parameter physical systems 
performing interrelated activities, at least one of such 

25 activities being designated as a risk- related activity 
and having an outcome level which may fall outside of 
boundary limits, said method developing a course of 
action (strategy) for said physical systems preventing 
any of said outcome levels for risk- related activities 

30 from falling outside of said boundary limits and 
comprising the steps of; 

(1) storing an initial model (11-17) of said 
systems in computer memory, said model including values 
for parameters of each of said physical systems; 

35 (2) producing a solution (19) of said initial 

model in said con^uter memory; 
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(3) deriving £rom said solution of said initial 
model a set of outcome levels (21-23) for said risk* 
related activities; 

(4) coznparing said set of outcome levels <25* 
5 29) for said risk- related activities with said boundary 

limits; 

(5) if the outcome of at, least one said risk* 
related activity falls outside said boundary limits for 
that activity, then forming a set of risk- limiting 

10 derivative constraints for all of said boundary limits; 

(6) esqpanding said initial model in said 
<:oznputer memory by adding said set of risk- limiting 
derivative constraints (31) to said initial model in said 
computer memory; 

15 (7) producing a solution (33) of said e3cpanded 

model with said set of added risk- limiting derivative 
constraints in said computer memory; 

(8) deriving from said solution (23, 35) of 
said esqpanded model a set of outcome levels for said 

20 risk-related activities; 

(9) comparing said set of outcome levels (25- 
29) for said risk- related activities with said boundary 
limits; 

(10) changing said esqpanded model by changing 
25 some of said boundary limits (37) and changing said risk- 
limiting derivative constraints (39) for said changed 
boundary limits; 

(11) producing a solution (33) of said changed 
model with said changed risk- limiting derivative 

30 constraints in said conputer memory; 

(12) deriving from said solution of said 
changed model a set of outcome levels (23, 35) for said 
risk - related activities ; 

(13) repeating steps (9) through (12) (23-25, 
35 29-39) xintil no outcome level of any of said risk-related 

activities falls outside said boundary limits for that 
activity; and 
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(14) Identifying a strategy (27) for said 
physical systems from the values of said strategic 
variables in said solution of said changed model. 

9. A con^uter method for managing risk in multiple 
5 parameter physical systems performing interrelated 
activities, at least one of such activities being 
designated as a risk- related activity and having an 
outcome level which may fall outside of boundary limits, 
said method developing a course of action (strategy) for 
10 said physical systems preventing any of said outcome 
levels for said risk- related activities from falling 
outside of said boundary limits and conprising the steps 
of: 

(1) storing an initial model (11-17) of said 
15 systems in said conputer memory, said model including 

both a set of decision variables and values for 
parameters of each of said physical systems, where 
strategic variables are a subset of said decision 
variables and where a strategy consists of a set of 
20 values specified for all said strategic variables; 

(2) producing a solution (19) of said initial 
model in said ccwnputer memory; 

(3) deriving from said solution of said initial 
model a set of outcome levels <23) for said risk-^^elated 

25 activities; 

(4) conparing said set of outcome levels (25- 
29) for said risk- related activities with said boundary 
limits and, if the outcome level of at least one said 
risk- related activity falls outside said boundary limits 

30 for that activity, then deriving at least one risk- 
limiting derivative constraint from such comparison; 

(5) expanding said initial model (31) in said 
computer memory by adding said risk- limiting derivative 
constraints to said initial model in said computer 

35 memory; 

(6) producing a solution (33) of said esqpanded 
model with said added risk- limiting derivative 
constraints in said computer memory; 
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(7) deriving from said soluCion (23, 35) of said 
expanded model a new set of outcome levels for said risk- 
related activities; 

<8) deriving from said solution of said 
5 expanded model a discretionary boundary limit (31) for at 
least one said risic- related activity; 

(9) comparing said new set of outcome levels 
(25-29) with said set of discretionary boundary limits 
and deriving at least one additional discretionary risk- 

10 limiting constraint from such comparison; 

(10) further expanding said expanded model (31, 
17-39) In said computer memory by adding said 
discretionary risk- limiting constraints to said expanded 
model In said conputer memory; 

15 (11) producing a solution (33) of said further 

expanded model with said discretionary risk- limiting 
constraints in said conqputer memory; 

(12) repeating steps (7) through (11) (23-25, 
29-39) until said solution of said further expanded model 

20 is satisfactory to decision -makers; and 

(13) identifying a strategy (27) for said 
physical systems from the values of said strategic 
variables in said solution of said expanded model. 

10. A computer method for managing risk in multiple 
25 parameter physical systems performing interrelated 
activities, at least one of such activities being 
designated as a risk -related activity and having an 
outcome level which may fall outside of boundary limits, 
said method developing a course of action (strategy) for 
30 said physical systeons preventing any of said outcome 
levels for said risk- related activities from falling 
outside of said boundary limits and comprising the steps 
of: 

(1) storing an initial symbolic mathematical 
35 model (11) of said systems in said coinputer memory, said 
model including a set of decision variables, where 
strategic variables are a subset of said decision 
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variables and where a strategy consists of a set of 
values specified for said strat glc variables; 

(2) entering values for parameters (13-15) of 
each of said physical systems Into said coi^puter memory; 
5 (3) converting said Initial symbolic 

mathematical model into an Initial computational model 
(17) in said computer memory by replacing symbols with 
numbers; 

(4) producing a solution (19> of said initial 
10 model in said computer memory; 

(5) deriving from said solution of said initial 
model a set of outcome levels (21-23) for said risk* 
related activities; 

(6) comparing said set of outcome levels (25) 
15 for said rleJc* related activities with said botmdary 

limits and, if the outcome level of at least one said 
risk-related activity falls outside said boundary limits 
for that activity, then deriving at least one rlek- 
limitlng derivative constraint from such conparlson; 
20 (7) expauuiing said initial model in said 

cOBiputer memory by adding (31, 37-39) said risk-ldlmiting 
derivative constraints to said initial model in said 
computer memory; 

(8) producing a solution (33) of said esqpanded 
25 model with said added risk-limiting derivative 

constraints in said cooputer memory; and 

(9) identifying a strategy (27) for said 
physical systems from the values of said strategic 
variables in said solution of said e^qpanded model * 

30 11. A computer method for managing risk under 

uncertainty in multiple parameter physi<»l systems 
performing interrelated activities, at least one of siwh 
activities being designated as a risk- related activity 
and having outcome levels which may fall outside of 

35 boundary limits, said systems having parameters with 

certain values and at least one parameter with uncertain 
values, said method developing a course of action 
(strategy) for said physical systems preventing any of 
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said outcome levels for said riflJc- related activities from 
falling outside of said boundary limits and comprising 
the steps of: 

(1) storing an initial model {115-119) of said 
5 systems in conputer memory, said model including both a 

set of decision variables and values tor parameters of 
each of said physical systems, where strategic variables 
are a subset of said decision variables and where a 
strategy consists of a set of values specified for said 
10 strategic variables, at least one of said parameters 

having a plurality of different possible values because 
of uncertainty I 

(2) producing solutions (121-127) of said 
initial model in said con^ter memory for a plurality of 

IS respectively different combinations of said parameters, 
where each of said different combinations includes a 
respectively different one of said possible values; 

(3) deriving from said solutions of said 
initial model a set of outcome levels (133) for said 

20 risk- related activities for each of said different 
combinations; 

(4) coi)q>aring said set of outcome levels (135) 
for said risk- related activities with said boundary 
limits and, if any of said boundary limits are violated, 

25 then deriving at least one risk- limiting derivative 
constraint from such comparison; 

(5) expanding said model (143) in said computer 
memory by adding said risk- limiting derivative 
constraints to said model in said ^onputer memory; 

30 (6) producing solutions (145) of said es^amded 

model with said added risk- limiting derivative 

constraints in said conqputer memory; 

(7) repeating iteratively steps (3) through (€) 

(129-137, 141-147) until solutions for all said different 
35 combinations of possible values are satisfactory for 

decision-makers; 
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(8) identifying a strategy (139) tor said 
physical systems from the values of said strategic 
variables in said solution of said expanded model ; and 

(9) recording outcome levels (149) in said 

5 ccsnputer memory for all scenario versus strategy versus 
ri8)c* related activity combinations into an outcome 
matrix. 

12. The computer risk management method of claim 11 
which includes the additional step of : 

10 (10) deriving from entries in said outcome 

matrix all scenario versus strategy versus risk- related 
activity combinations of regret values and forming a 
regret matrix (FIG. 8) in said conq>uter memory, where 
regret values, for each risk* related activity, are 

15 computed by the equation 

REGRET (H](J1 « BEST (HI - 0(H) (JJ , 
where BBST[H] is the best scenario outcome for scenario H 
and 0[H1 fJ] is the outcome for strategy J under scenario 

20 13. The computer risk management method of claim 12 

in which each said candidate strategy is described in a 
multidimensional format by the outcome and regret values 
for said risk- related activities, such 
multidimensionality taking into consideration factors 

25 including the multiplicity of scenarios, risk- related 

activities, and causes of uncertainty, and which aids in 
selection of a strategy on the basis of such 
multidimensional description. 

14. The conqputer risk management method of claim 12 
30 in which both of said matrices have n dimensions, where n 

is an integer greater tham one. 

15. A computer method for managing risk under 
uncertainty in multiple parameter physical systems 
performing interrelated activities, at least one of such 

35 activities being designated as a risk- related activity 
and having an outcome level which may fall outside of 
boundary limits, said systems having parameters with 
certain values and at least one parameter with uncei?tain 
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values, said method developing a course o£ action 
(strategy) for said physical systems preventing auiy of 
said outcome levels for said risk- related activities from 
falling outside of said boundary limits and comprising 
5 the steps of: 

(1) storing an initial model (115- 119) of said 
systems in conputer memory « said initial model including 
values for parameters of each of said physi^cal systems, 
at least one of said parameters being uncertain in that 

10 it has a plurality of different possible values because 
of uncertainty, and said model comprising interconnected 
sukmiodels for a plurality of respectively different 
coaibinatlons of said parameters, where each of said 
different combinations is a respectively different 

15 scenario which Includes a respectively different one of 
said possible values; 

(2) entering values for ^ach of said certain 
parameters (119) into all of said scenario submodels in 
said cooputer memory; 

20 (3) entering a respectively different one of 

said possible values (121-123) for said uncertain 
parameters into each of said scenario submodels in said 
cc»nputer memory; 

(4) producing a solution (127) of said model in 
25 said computer memory, said solution including a 

subsolution for each of said scenario submodels; 

(5) deriving from said subsolutlons for each of 
said scenario submodels respective outcome levels (133) 
for said ris)c- related activities f or -^ach said scenario; 

30 (6) conaparing said set of outcome levels (135) 

for said risk- related activities from each of said 
scenario submodels with said boundary limits and, if any 
of said boundary limits are violated, then deriving at 
least one risk- limiting derivative constraint from such 

35 conparison; 

(7) e3q;>anding at least some of said scenario 
submodels (143) in said computer memory by adding said 
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risk- limiting derivative constraints to such expanded 
submodels; 

(8) producing solutions of said expanded 
submodels (145) with said added risk-limiting derivative 

5 constraints in said conqputer memory; 

(9) repeating iteratively steps (5) through (8) 
(133-137, 141-147) until said solutions for all said 
different scenarios are satisfactory for decision-makers; 
and 

10 (10) identifying a strategy (131) -for said 

physical systems from the values of said strategic 
variables in said solution of said expanded model. 

16. The cooputer risk management method of claim 15 
in which said discretionary limits are imposed on said 

IS outcome and regret levels of said risk-related activit;ies 
under specified individual scenarios or soenario groups 
in said solutions of said computational model. 

17. A computer method for managing risk under 
uncertainty in multiple parameter physi<:al systems 

20 performing interrelated activities, at least one of such 
activities being designated as a risk- related activity 
and having an outcome level which may fall outside of 
boundary limits, said systems having parameters with 
certain values and at least one parameter with uncertain 

25 values, said method developing a course of action 

(strategy) for said physical systems preventing any of 
said outcome levels for said risk-related activities from 
falling outside of said boundary limits, said method 
comprising the steps of: 

30 (1) storing an initial model (41-45) of said 

physical systems in computer memory, said imodel including 
both a set of decision variables and values for 
parameters of each of said physical systems, at least^ one 
of said paurameters being uncertain in that it has a 

35 plurality of different possible values because of 

uncertainty, and said model comprising single -scenario 
computational models (49-53) f or jeach respectfully 
different combination of said parameters and each of said 
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different conibinations is a scenario, said single- 
scenario cosiputational models including constraints 
arising from predetermined boundary limits, and strategic 
variables are similar subsets of decision variables in 
5 each of said single -scenario coiniutational models, where 
the possible range of values for each said strategic 
variable is divided into segments (47) , and where a 
strategy consists of a set of segments tor specified 
values of said strategic variables; 

(2) producing solutions <55) of said single- 
scenario computational models in said computer memory €or 
all scenarios, deriving from said solutions of said 
initial model a set of said outcome levels (57) tor said 
risk- related activities for each of said scenarios, and 

15 comparing said outcome levels (59) for said risk-related 
activities with said boundary limits; 

(3) if at least one of said scenario solutions 
is infeasible, in that at least one outcome falls outside 
predetermined boundary limits for a risk- related 

20 activity, then modifying said boundary limits (67) until 
all said scenario solutions become feasible; 

(4) additionally imposing discretionary 
constraints on outcome levels (€3) of said risk- related 
activities in specified single- scenario coii5)utational 

25 models and modifying said discretionary constraints until 
all of them are satisfied under all said scenarios; 

(5) analyzing the satisfactory solutions (75) 
of said single -scenario computational models for all 
scenarios, finding the segments into which the values of 

30 each strategic variable in each said scenario solution 

fall (77), and forming a candidate strategy (79) for each 
subset of said solutions that has the same combination of 
strategic variables in specific segments over the whole 
set of said strategic variables, while subsets of 

35 solutions with respectively different combinations form 
different candidate strategies; 

(6) forming strategy definition constraints 
(87) for a candidate strategy, where said strategy 
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definition constraints ensure that the value of ^ach 
strategic variable falls into a specified segment, adding 
said strategy definition constraints (91) in said 
contputer memory in turn to said single -scenario 
5 con^utatlonal models for all of said scenarios, and 

generating solutions (93) in said computer memory for all 
of said expanded single -scenario computational models; 

(7) if at least one solution of said single - 
scenario contqputatlonal models is infeeuslble, -^hen 

10 excluding said csuididate strategy (97) fxom further 
cons iderat ion ; 

(8) If said candidate strategy remains viable, 
then recording outcome levels (109-111) for all scenario 
versus strategy versus rls)c*related activi-ty combinations 

15 into an outcome matrix; and 

(9) repeating steps (6) through (B) (87-97, 
109-111) for the whole set of candidate strategies formed 
in step (5) . 

18. The ccmqputer method of claim 17 in whi^ a 

20 plurality of said scenarios are grouped (161) before 
forming single- scenario computational models in said 
computer memory and in which single -scenario 
computational models In said computer memory are formed 
only for said grouped scenarios (165) . 

25 19 . The '<:on9>uter method of claim 18 In which a 

respectively different plurality of scenarios and single- 
scenario computational models formed therefrom are 
regrouped (161) iteratively In said computer memory. 
20, A computer method for managing risk under 

30 uncertainty in multiple parameter physical systems 

performing interrelated activities, at least one of such 
activities being designated as a risk-«lated activity 
and havdlng an outcome level which may fall outside of 
boundary limits, said systems having parameters with 

35 certain values and at least one parameter with uncertain 
values, said method developing a ^course of action 
(strategy) for said physical systems preventing any of 
Bald outcome levels for said risk- related activities from 
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falling outside of said boundary limits and conqprising 
the steps of: 

(X) storing an initial model (llS-123) of said 
physical systems in contputer memory, said model including 
5 both a set of decision variables and values for 

parameters of each of said physical systems, at least one 
of said parameters being uncertain in t:hat it has a 
plurality of different possible values because of 
uncertainty, said model being a multiscenario 

10 contputational model that conprlses interconnected 

submodels for each respectively different combination of 
said parameters and each of said different combinations 
is a scenario, said scenario submodels including 
constraints arising from said predetermined boundary 

15 limits and strategic variables are similar subsets for 

each scenario of said decision variables of said scenario 
submodels, where the possible range of values for each 
said strategic variable is divided into segments, where a 
strategy consists of a set of segments for specified 

20 values of said strategic variables, and where under each 
scenario the values of each strategic variable €all into 
the same segment; 

(2) producing a solution (127) of said 
multiscenario conqputational model in said coiiq>uter 

25 memory, said solution conqprising scenario subsolutions 

for all scenario submodels, deriving from said solutions 
of said submodels a set of outcome levels for said risk- 
related activities for each of said scenarios, and 
con^aring said outcome levels for said rislc-related 

30 activities with said boundary limits; 

(3) if at least one of said scenario 
subsolutions is infeasible in that at least one outcome 
falls outside said predetermined boundary limits for a 
risk-related activity in said scenario, then modifying 

35 said boundary limits until all of said scenario 
subsolutions become feasible; 

(4) additionally is^osing discretionary 
constraints (143) on outcome levels of said risk-related 
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activities in specified scenario siifamodels and inodi£ying 
said discretionary constraints until all of them are 
satisfied in all of said scenario siibmodels; 

(5) analyzing satisfactory solutions of said 
5 scenario submodels for all scenarios, finding segments 
into which the values of each strategic variable fall, 
and forming a candidate strategy (131) that reflects the 
combination of strategic variables in specific segments 
obtained over the whole set of strategic variables; 
10 <6) recording outcome levels (133, 149) tor all 

scenario versus strategy versus risk- related activity 
conibinations into an outcome matrix; 

(7) modifying said discretionary constraints on 
outcome levels (143) for specified risk-^^elated 

15 activities and scenarios to obtain another feasible and 
satisfactory candidate strategy; and 

(8) repeating steps (5) chrough (7) (131-137, 
143-149) until a sufficient number of -candidate 
strategies is obtained. 

20 21. The computer method of claim 20 in whi<:h a 

plurality of said scenarios are grouped (161) before 
forming scenario submodels in said multiscenario 
con^utational model in said computer memory cuid in which 
scenario submodels in said cosputer memory are formed 

25 only for said grouped scenarios. 

22. The computer method of claim 21 in which a 
respectively different plural it:y of said scenarios and 
scenario submodels formed therefrom are regrouped (161) 
iteratively in said cooqputer memory* 

30 23. A computer method for managing risk under 

uncertainty in multiple parameter physical systems 
performing interrelated activities, at least one of such 
activicies being designated as a risk- related activity 
and having an outcome level which may fall outside of 

35 boiuidary limits, said systmis having parameters with 

certain values and at least one parameter with uncertain 
values, said method developing a course of action 
(strategy) for said physical syst ms preventing any of 
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said outccsme levels for said risk -related activities from 
falling outside of said boundary lixnlts and comprising 
the steps of: 

(X) storing an initial model (191-195) of said 
5 physical systems in coznputer memory, said model including 
both a set of decision variables and values for 
parameters of each of said physical systems « at least one 
of said parameters being uncer^iain in that it has a 
plurality of different possible values because of 

10 uncertainty, and said model being a con?>utational model 
that comprises multiple group multiscenario models for a 
multitude of groups of each respectively different 
combination of said parameters, where each of said 
different coxnbinations is a scenario and each of said 

15 multiple groups cosprises a plurality of said scenarios, 
said group multiscenario models including constraints 
arising from predetermined boundary limits, and strategic 
variables are similar subsets of decision variables for 
each of said group multiscenario models, where the 

20 possible range of values for each strategic variable is 

divided into segments, and %^ere a strategy consists of a 
set of segments for specified values of said strategic 
variables; 

(2) grouping respective pluralities of 
25 scenarios into scenario groups (199) i 

(3) forming a respective group multiscenario 
computational model (205) for each scenario group; 

(4) producing solutions of said group 
multiscenario computational models (207) in said conqpiiter 

30 memory, said solutions comprising scenario subsolutions 
for all scenario submodels, deriving from said solutions 
of said submodels a set of said outcome levels (209) for 
said risJc- related activities for each of said scenarios, 
and conparing said outcome levels (211) for said risk- 

35 related activities %ri.th said boundary limits; 

(5) if at least one of said scenario 
subsolutions is inf easible in that at Least one outcome 
falls outside said predetermined boundary limits for a 
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rlak- related activity in its scenario * then modifying 
(221) said boundary limits until all said scenario 
subsolutions become feasible; 

(6) additionally litposing discretionary 

5 constraints (219) on outcome levels of said risk-related 
activities in specified scenario sxibmodels and modifying 
said discretionary constraints until all of them are 
satisfied in all of said scenario submodels; 

(7) repeating steps (3) through (6) (20S-222) 
10 for the whole set of scenario groups; 

(8) analyzing said satisfactory solutions of 
said scenario submodels for all said scenarios, finding 
which of said scenarios have worst outcomes (223). 
forming a new scenario group with the worst outcomes, 

15 forming a new group multiscenarlo contputational model 
(224) for that group and producing a solution (225) of 
said new multiscenarlo computatiozial model; 

(9) if at least one of said scenario 
subsolutions in said new model is Infeasible. then repeat 

20 steps (5) and (6) (219-221) for the new group 
multiscenarlo computational model; 

(10) analyzing into which segments said values 
of each strategic variable fall (231). and forming a 
candidate strategy (239) that reflects the obtained 

25 combination of strategic variable in specific segments 
over the whole set of said strategic variables; 

(11) forming strategy definition constraints 
(239) for said candidate strategy; 

(12) adding said strategy definition 

30 constraints to the group multiscenarlo computational 

model (243) for the initial scenario group and solving 
the resulting group multiscenarlo computational model; 

(13) if at least one of said scenario 
stibsolutions Is infeasible, then screening out said 

35 candidate strategy (247, 2S1-25S) , modifying the 

coR^sltlon of the worst scenario group, and returning to 
step (8) ; 
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(14) recording out:come levels (257-259) for all 
scenario versus strategy versus rislc- related activity 
coxnbinations for all scenarios of said scenario group 
into an outcome matrix; 
5 (15) repeating steps (11) through (14) (239- 

259) for the whole set of said scenario groups; and 

(16) repeating steps (8) through (15) (223-259) 
until a plurality of candidate strategies are obtained. 
24. The computer method of claim 23 in whi^h a 

10 plurality of said scenarios are grouped (1^1) before 
forming scenario submodels in said multiscenario 
computational models in said computer memory anA in which 
scenario sulnnodels in said computer memory are formed 
only for said grouped scenarios. 

15 25. The computer method of claim 24 in which a 

respectively different plurality of said scenarios and 
scenario submodels formed therefrom are regrouped (161) 
iteratively in said con^uter memory. 

26. A computer method for managing rislc under 

20 uncertainty in multiple parameter physical systems 

performing interrelated activities, at least one of such 
activities being designated as a ris)c- related activity 
and having an outcome level which may fall outside of 
boundary limits, said systems having parameters with 

25 certain values and at least one parameter with uncertain 
values « said method developing a course of action 
(strategy) for said physical systems preventing any of 
said outcome levels for said risk- related activities from 
falling outside of said boundary limits and conprlsing 

30 the steps of: 

(1) storing an initial model (191-195) of said 
physical systems in computer memory, said model including 
both a set of decision variables and values for 
parameters of each of said physical systems, at least one 

35 of said parameters being uncertain in that it has a 
plurality of different possible valxies because of 
uncertainty, and said model being a coxnputational model 
that conqprises multiple group multiscenario models for a 
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multitixde of groups o£ each respectively different 
combication of said parameters « where each of said 
different combinations is a scenario and each of said 
multiple groups comprises a plurality of said scenarios « 
5 said group nultiscenario models including constraints 

arising from predetermined boundary limits, and strategic 
variables are similar subsets of decision variables for 
each of said group multiscenario models, where the 
possible range of values for each strategic variable is 
10 divided into segments, and where a strategy consists of a 
set of segments for specified values of said strategic 
variables; 

(2) grouping respective pluralities of 
scenarios into scenario groups (199) ; 

15 (3) forming a respective multiscenario 

CQn^utational model (205) for each scenario group; 

(4) producing solutions of said multiscenario 
computational models (207) in said -conqputer memory, said 
solutions coziqprising scenario subsolutions for all 

20 scenario submodels, deriving from said solutions of said 
submodels a set of said outcome levels (209) for said 
risic- related activities for each of said scenarios, and 
conparing said outcome levels (211) for said risk- related 
activities with said boundary limits; 

25 (5) if at least one of said scenario 

subsolutions is Inf easlble in that at least one outcome 
falls outside said predetermined boundary limits for a 
risk- related activity in its scenario, then modifying 
(221) said boundary limits until all said scenario 

30 subsolutions become feasible; 

(6) additionally imposing discretionary 
constraints (219) on outcome levels of said risk- related 
activities in specified scenario submodels and modifying 
said discretionary constraints until all of them are 

35 satisfied in all of said scenario sxibmodels; 

(7) repeating steps (3) through (6) (20S-232) 
for the whole set of scenario groups; 
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(8) analyzing said .satisfaccoxy 'solutions of 
said scenario submodels £or all said scenarios, finding 
which of said scenarios have worst out<:ome8 (233), 
forming a new scenario group with the worst outcomes, 

5 forming a nev multiscenario computational model (234) for 
that group and producing a solution (22S) of said new 
multiscenario computational model; 

(9) if at least one of said scenario 
Bubsolutions is infeasible, then repeat steps (5) and (€) 

10 (219-221) for the new multiscenario coinputational model; 

(10) analyzing into which segments said values 
of each strategic variable fall (231), and forming a 
candidate strategy (239) that reflects the obtained 
conibination of strategic variable in specific segments 

15 over the whole set of said strategic varial>les; 

(11) forming strategy definition constraints 
(239) for said candidate strategy; 

(12) adding said strategy definition 
constraints (243) to the multiscenario computational 

20 model for the initial scenario group and solving the 
resulting multiscenario computational model; 

(13) if at least one of said scenario 
subsolutions is infeasible, then screening out said 
candidate strategy (247, 251-255), modifying the 

25 composition of the worst scenario group, and returning to 
step (8) ; 

(14) recording outcome levels <257-2S9) for all 
scenario versus strategy versus rlslc- related activity 
<:Qmbinations for all scenarios of said scenario group 

30 into an outcome matrix; 

(15) repeating steps (11) through (14) (239- 
259) for the whole set of said scenario groups, where 
tightening of said strategy definition constraints is 
derived from solutions of said group multiscenario models 

35 for worst scenario groups and is used as a decostposition 
method to solve the large underlying overall 
multiscenario model; and 
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(16) repeating steps (8) through (15) (223-259) 
until a plurality of candidate strategies are obtained. 

27. The ccKRputer rislc management method of claim 12 
in which the simultaneous application of a multitude of 
5 conversion weights and coefficients to combine and 

converse entries of out;coroe and regret, matri-ces is used 
(269-279« 285-295) as an approximation technique to deal 
with fflultidimensionality of matrices « multiple optimality 
criteria, and qualitative outcomes. 
10 28. The conqputer risk management method of claim 6 

in which comparison of said candidate strategies is based 
on three jointly used regret -based synthetic criteria for 
two-dimensional outcome (payoff) and regret matrioea: 

(a) OPR * The Optimism- Pessimism Index (Regret) 
15 Criterion (PIG. 9); 

(b) PIR - The Partial ignorance (Regret) Criterion 
(FIG. ID) ; and 

(c) MIRR * The Modified Insufficient Reason (Regret) 
Criterion (FIG. 11) / where 

20 said OPR criterion is: 

Select J that minimizes 
ZOPR(J,AIiPHIl) " (1 - ALPHA) ♦ BEST REGRET (J) ♦ 
ALPHA * WORST RBGRBT<J) , 
0 ;i ALPHA ft 1 ( J»l, . . . , NST) , (1) 

25 the value of ALPHA is specified by the user and NST is 
the numiber of strategies; BEST REGRET (Or) and WORST 
REGRET (J) are the best and worst values of regret for 
strategy J; and regret in Eq. (1) is coxnputed by ^q. i2) : 
REGRET (HI I J] « BESTtH] - PAYOFF [H] IJ] , (2) 
30 where BEST(H) is the best payoff for any strategy under 
scenario H and PAYOFF [H] [J) is the payoff for stratiegy J 
under scenario H; 

said PIR cricerion is: 

Select J that minimizes 
35 ZPIR(J,ALPHA) » (1 - ALPHA) ♦ EXPECTOD REGRET (J) ♦ 

ALPHA * WORSTs REGRET (J) , 

0 M ALPHA II 1 (J-1,..., NST)^ (3) 
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the value of ALPHA is again specified by the user, and 
EXPECTED REGRET (J) i8 computed by Bq. (4): 

S 

EXPECTED RBGRBT (J) - 0 u IHl REGRET (H) ( JJ , (4) 

5 

REGRET [H] [J] is regret for strategy J under scenario H, 
u{H] is the probability of scenario H, and S is the 
nuniber of scenarios; and 

said MIRR criterion is: 

Select J that minimizes 

ZIRR{J,ALPHA) - (1 - ALPHA) ♦ EXPECTED IR REGRET (J) 

ALPHA * WORST REGRET (J) , 
0 M ALPHA fill 1 {J«l, . • . #NST) , <5) 

15 the value of ALPHA is again specified by the user, and 
EXPECTED IR REGRET (J) is computed by Eq. (6) : 

S 

EXPECTED IR REGRET (OT) « REGRET [HI £J1 > /S- (6) 

H«l 

20 29. The can?juter risk management method of -claim 28 

in which coa5>arison of said strategies is based on three 
jointly used regret -based 'efficient frontier" methods 
(FIGS. 8, 9, 10) which are respectively based on said 
CPR, PIR and MIRR synthetic criteria for two-dimensional 

25 regret matrices; said efficient frontier is the set of 

best strategies on the [0, U range of probability ALPHA; 
said OPR efficient frontier cari5>rises straight- line 
segments defined by Eq. (l) , said PiR efficient frontier 
comprises straight-line segments defined by Eq. (3), and 

30 said MIRR efficient frontier con^rises straight- line 
segments defined by Eq. (5), respectively; 

where each strategy on the .efficient frontier 
is characterized by a subinterval of the [0, 1] range ^n 
which this strategy is the c^timal choice, and the 

35 efficient frontier incorporates significant information 
about the relative merits and demerits of any strategy: 
• the composition of the subset of strategies that 
form the efficient frontier; 
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• the width of the subinterval supporting each 
optixnal strategy; 

• the order o£ optimal strategies from the 
optimistic end of the [0, 11 interval to the pessimistic 

5 end; and 

• the difference between the regrets for the best 
strategy and other strategies, which shows the possible 
impairment of results in choosing a non- optimal strategy. 

30. The conputer risk management method of claim 29 
10 in which the three OPR, PIR and MIRR synthetic criteria 

(FIGS. B, 9, 10) for said two-dimensional regret matrices 
and the three regret -based OPR, PIR and MIRR "efficient 
frontier" methods are jointly used to find desirable 
limits on the tightening of discretionary boundaries and 
15 associated discretionary constraints . 
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